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ABSTRACT 

The XMM-Newton - Blanco Cosmology Survey project (XMM-BCS) is a coordinated X-ray, optical and mid-infrared cluster survey in 
a field also covered by Sunyaev-Zel'dovich effect (SZE) surveys by the South Pole Telescope and the Atacama Cosmology Telescope. 
The aim of the project is to study the cluster population in a 14 deg^ field (center: a « 23:29:18.4, 6 ~ -54:40:33.6). The uniform 
multi-wavelength coverage will also allow us for the first time to comprehensively compare the selection function of the different 
cluster detection approaches in a single test field and perform a cross-calibration of cluster scaling relations. 

In this work, we present a catalog of 46 X-ray selected clusters from the initial 6 deg^ survey core. We describe the XMM-BCS source 
detection pipeline and derive physical properties of the clusters. We provide photometric redshift estimates derived from the BCS 
imaging data and spectroscopic redshift measurements for a low redshift subset of the clusters. The photometric redshift estimates are 
found to be unbiased and in good agreement with the spectroscopic values. 

Our multi- wavelength approach gives us a comprehensive look at the cluster and group population up to redshifts z ~ I- The median 
redshift of the sample is 0.47 and the median mass Msoo ~ 1 x 10'* Mq (~ 2 keV). From the sample, we derive the cluster log A' -log 5 
using an approximation to the survey selection function and find it in good agreement with previous studies. 

We compare optical mass estimates from the Southern Cosmology Survey available for part of our cluster sample with our estimates 
derived from the X-ray luminosity. Weak tensing masses available for a subset of the cluster sample are in agreement with our 
estimates. Optical masses based on cluster richness and total optical luminosity are found to be significantly higher than the X-ray 
values. 

The present results illustrate the excellent potential of medium-deep. X-ray surveys to deliver cluster samples for cosmological 
modelling. In combination with available multi-wavelength data in optical, near-infrared and SZE, this will allow us to probe the 
dependence of the selection functions on relevant cluster observables and provide thus an important input for upcoming large-area 
multi-wavelength cluster surveys. 

Key words. Surveys, Catalogs, Galaxies: clusters: general. Cosmology: Large-scale structure of Universe 



1. Introduction 

The formation of the cold dark matter (CDM) dominated large- 
scale structure of the Universe is hierarchical with smallest ob- 
jects collapsing first. With passing time more and more massive 
structures are able to decouple from the Hubble flow and en- 
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ter the non-linear regime, collapse and eventually virialize. The 
statistical properties of the matter density field (e.g. its power 
spectrum) as well as the growth of the structures are strongly de- 
pendent on the background cosmology and can be thus used to 
put constraints on cosmological models. 

From this point of view, clusters occupy a very important 
place in the structure formation scenario, by being the most re- 
cent (i.e. redshifts z < 2 - coincident with the onset of the dark 
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energy dominance) and thus also the most massive structures 
(lO'"^ - lO'^ M0) to virialize. The cluster abundance is therefore 
exponentially sensitive to the growth of the large s cale-structure 
and to the underlying cosmological parameters (|Haiman et alJ 
l2001l;lMaiumdar & Mohiil2003l:lHaiman et alJl2005lK 

The key parameter in cosmological tests of this type - the 
total mass of clusters (identified with dark matter halos) - is it- 
self not directly observable. Fortunately, in first approximation, 
clusters are virialized and their growth is gravitationally driven 
and therefore self-similar This allows us to link their mass to 
some suitable observable quantity originating from the baryonic 
components of a cluster - its galaxy population and the intra- 
cluster medium (ICM). The ICM is directly observable in X-rays 
or through the distortion of the Cosmic Microwave Background 
(CMB) imprinted by the ICM thermal electron population via 
inverse Com pton scattering (the so- called Sunyaev-Zel'dovich 
effect (SZEl lSunvaev & Zerdovic"hl[T97l . 

Since the ICM closely traces the DM potential, it offers 
better (i.e. lower scatter) mass-proxies than those available 
from optical observati ons of the cluster's galaxy population 
(e.g. iReves et al.l r2008h . In X-rays, the simplest and observa- 
tion ally least expensive mass-p r oxy is the X-ray luminosity 
Ly (iReiprich & BohringeJ 120021 iPratt et al.ll2009t iMantz et al.1 
I2010ar 

For the SZE experiments the most direct way to esti- 
mate the cluster mass is from the source signal-to-noise ra- 
tio (e.g. Williamson et al. 2011; Vanderlinde et al. 2010) and 
more importantly, through the integrated Compton parameter 
Ysz- Numerical simulat ions suggest that F s z is an excellent 
proxy of cl uster mass (Ida Silva et alj 120041: iMotl et al.1 120051; 
lNagaill2006l) . First cross-comparisons with X-ray and SZE stud- 
ies are generally finding good agreement between the mass 
estimates a nd no significant deviation fro m the self-similar 
predictions (Planck Collaboratio nll201 Idlbllal; iMelin et a l.ll201 U 
Andersson et al. 2010; Marron ^et al.l 20091 iBonament e et aP 
20081) . 



If deeper X-ray observations are available, we can use the 
spectroscopic temperature Tx, gas mass Mg and their combi- 
nation Yx = TxM„ (t he X-ray a nalogue t o the Fsz parameter , 
iKravtsov et a l. 200 6: Vikhlinin et al. 20091 lArnaud et al.ll2010l) 
as goo d mass proxies. Using the Yx parameter I Vikhhnin et al.l 
( I2OO9 *) put a strong constraint on the cosmological parameters 
including the dark energy equation of state. From a methodolog- 
ical point of view, this is interesting for two reasons: 1) it shows 
that useful cosmological constrains can be obtained already from 
relatively small samples of clusters of galaxies, demonstrating 
the exceptional potential of this type of cosmological tests; and 
2) already this modest sample is practically systematics-limited, 
especially due to uncertainties in the mass estimation. 

There are many factors that affect the scaling relations 
and the intrinsic scatter of the cluster populations around 
these relations: the presence of cool cores (Markevitch 1998; 
lO'Hara et al.ll2006tlMotl etal.l2005; Prat t et al. 200 9), substruc- 
tures and the cluste r's dynamical state feohringe r et al.ll20Tot 
iJeltema et al.' l2008l) and additional non-gravitational physics 
(Nagai 2006), etc. In addition, one has to account for the 
Malmquist and Eddington bias when determining the scaling 
relations from an X-ray selected sample of clusters by proper 
trea tment of the select i on and mass functions (especially for 
Lx, 'Pacaud et al."2007t I Vikhhnin et aLll2009t iPratt et al.ll2009t 
[Mantz et al. 2010a b). As our cluster samples cover broader red- 
shift ranges potential deviations from the self-similar evolution 
of the scaling relations also become an important question. 



In summary, in order to be able to well constrain cosmolog- 
ical models with cluster samples we need: 1) large cluster sam- 
ples covering redshifts beyond unity; 2) good knowledge of the 
cluster selection function's dependence on relevant observables 
and the distributions of these observables in the cluster popula- 
tion; 3) a reliable, low scatter mass-proxy with a known evolu- 
tion in the redshift range of interest. 

Surveying for clusters in SZE has a large potential with re- 
gards to all three requirements, having an almost redshift in- 
dependent selection very close to a selection function with a 
fixed mass limit at all redshifts and a robust mass-proxy in the 
Fsz parameter. Two ground-based large-area cluster surveys are 
currently underway: one by the South Pole Telescope (SPT) 
and one by the Atacama Cosmology Telescope (ACT). Both 
have already provided their first SZE-selected cluster s amples 
(Williamson et al. 201 1; Vanderlinde et al. 2010; Marria ge et al.l 
120 lot IStaniszewski et a l. 2009) as well a s observations of al - 
ready known clusters (iPlagge et 10120101: iHincks et alJl20Tol) . 
Also the Planck space mission ha s delivered its first cluster cat- 
alog jPlanck Collaborationl201 lah . 

While the SZE surveying approach is a very interesting new 
channel to perform cluster cosmology, there is still much work 
to be done at these early stages to understand the systematics 
like the influence of radio/sub-mm sources and primary CMB 
fluctuations on the selection, the mass calibration and sensitivity 
to cluster outskirts. 

A multi-wavelength follow-up program of SZE selected 
clusters is essential, but selection function studies require also 
comparison of blind surveys. To this end we are conducting the 
XMM-BCS cluster survey. The survey field covers a 14 deg^ 
area in the overlap region of the SPT and ACT surveys. The field 
has full coverage with the 4m CTIO telescope at Cerro Tololo, 
Chile, in the framework of the Blanco Cosmology Survey (BCS) 
in griz bands and Spitzer observations in the mid-infrared (mid- 
IR). With this optical to mid-IR coverage we are able to provide 
robust photometric redshift estimates out to redshifts ^ 0.8 (~ 1 
once also the Spitzer data is included). The X-ray coverage con- 
sists of XMM-Newton observations split into two distinct parts. 
The 6 deg^ core of the X-ray survey field was observed with 42 
individual, standard pointings (with ~ 10 ks effective exposure 
time). In this work, we present an initial cluster catalog based on 
these observations. 

After SPT commenced its operations, it was soon found that 
the mass threshold of contemporary SZE surveys is higher than 
expected. In order to offer a larger overlap between the SZE and 
X-ray selected cluster samples, we carried out an extension of 
the X-ray survey by covering an additional 8 deg^ in three large- 
area fields utilizing the new mosaic mode type of observations. 
These observations allowed us to cover a significantly larger area 
in a very time-efficient way. First results as well as details on the 
analysis of this type of XMM-A^ewfon observations are described 
in ISuhada et al.l (l2010l) . We demonstrate there the feasibility of 
blindly detecting clusters found with current generation SZE ex- 
periments in only ~ 3 ks long XMM-Newton observations (in- 
cluding tentative spectroscopic temperature measurements) in 
the case of two SPT detected clusters. The final 14 deg^ X-ray 
cluster catalog is expected to roughly double the number of clus- 
ters in the present sample and this sample will then be interesting 
also for its cosmology-constraining power 

The paper is organized as follows: in Sect. |2] and |3] we de- 
scribe the analysis of the X-ray observations and cluster detec- 
tion pipeline. The optical data, photometric redshift estimation 
and spectroscopic campaign are detailed in Sect. |4| In Sect. |5] 
we provide our cluster sample, the physical parameters of the 
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detected clusters and determine the survey's preliminary sta- 
tistical properties. We also cross-correlate our cluster catalog 
with known sources and carry ou t a detailed comparison with 
the optically sel e cted sample of iMenanteau et alJ (l2009l) and 
iMenanteau etal] (|201(I) (M09 and MIO hereafter). Sect. |6] dis- 
cusses the X-ray error budget and gives an outlook on the up- 
coming work in the context of the XMM-BCS survey. We give 
our conclusions in Sect. |5] In the appendices we provide ancil- 
lary information for the individual clusters, a preliminary com- 
parison of our simplified sensitivity function calculations with 
realistic simulations and a cross-comparison with the XMM- 
LSS cluster survey. 

We adopt a ACDM cosmology with (Oa, w, i/o) - 
(0.7, 0.3, -1, 70 km s"' Mpc"'). Estimated physical parameters 
are given in apertures corresponding to overdensities by factors 
200 and 500 with respect to the critical density of the Universe 
at the redshift of a given cluster Throughout the article we refer 
to objects in our sample as "clusters" regardless of their mass. 
The term "group" will be used to refer to systems with masses 
^200 S 10'^ We will refer to individual objects by their 
identification number (ID). Proper object names are listed in 
Table O 

2. XMM-A/eivton data reduction 

The XMM-Newton coverage of the XMM-BCS survey core con- 
sists of 42 partially overlapping pointings with offsets of 22.8 
arcmin covering a total area of about 6 deg^ (see Fig. [T]). The 
observing time was allocated in the frame of an XMM-Newton 
Large Program during A06. Four additional observations were 
carried out in A07 to replace fields with large losses due to soft- 
proton flaring. The observation of field F09 (Table [U was car- 
ried out in two parts. The total observing time amounts to ~ 580 
ks, with an average total nominal time per pointing of ~ 15 ks 
(including instrument setup time and high background periods). 
Table[T]displays the basic information about the individual point- 
ings. The THIN filter was used in all observations. The EPIC PN 
camera was operated in full frame mode. 

The full XMM-BCS X-ray field is displayed in Fig. [T] The 
core region presented in this work is inside the white boundaries 
(region F). Regions A, B and C mark the three mosaic extensions 
of the survey. The five missing fields in region F have been com- 
pletely lost due to flaring (F03, F05, F42) or had large time losses 
due to flaring and have a very high residual quiescent soft pro- 
ton contamination (F07 and FI3). The point source subtracted 
signal-to-noise map of the core region is displayed in Fig.|2](see 
e.g. Finoguenov et al. 2009; Bielbv et al. 2010). 

The EPIC data was processed with the XMM-Newton 
Standard Analysis System (S AS) version 7. 1 .CQ- We reduced and 
calibrated the raw observational data files with the SAS tasks 
epchain for the EPIC PN detector and emchain for both MOS 
detectors. Events in bad pixels, bad columns and close to the 
chip gaps are excluded from further analysis. 

The event lists were screened for high background periods 
cause d by soft proton flares with a two-step sigma clipping algo- 
rithm (Pratt & Arnaud 2003). We reject time intervals with back- 
ground count rates above the 3cr limit from the mean level in the 
12 - 14 keV band for PN and 10 - 12 keV band for MOSl and 
MOS2. The mean background count rate is determined by fitting 
a Gaussian model to the distribution of counts in the light curve 
binned in 100 second intervals. After this first cleaning step, we 
apply the same 3cr clipping procedure in the 0.3 - 10 keV band 

' We provide a test using the current SAS 11. 0.0 in Sect. 16. II 



on 10 second binned light curves to conservatively remove time 
intervals affected by low energy flares. An example of a two-step 
cleaned light curve is displayed in Fig. [3] 

Time lost due to flaring in our observations amounts typi- 
cally to ~ 20% of the full effective observing time. Six observa- 
tions of the initial fields from A06 were too heavily affected by 
the flaring even after the two step cleaning. Three of these fields 
have been replaced by observations in A07 (FOlb, F02b, F35b); 
the partially lost field F04 was reobserved as well. 

Detection and analysis of faint diffuse sources like clusters 
of galaxies in shallow surveys can be additionally affected by 
low energy soft protons with a roughly constant flux. This so- 
called quiescent soft proton background can not be detected 
based on light curve screening due to its small temporal vari- 
ations, especially not in observations with a short duration. In 
order to characterize possible contamination from this part of the 
non-X-ray background, we applied the diagnostics developed by 
iDe Luca & Molendil (l2004 . based on flux ratios inside and out- 
side the field of view of each detector. The vast majority of fields 
is not contaminated by the quiescent soft proton background at 
all in any of the detectors. Four fields (F04, F06, F16, F25) have a 
slighj3 contamination with negligible effect on data analysis and 
derived results. Fields F07, F13 have significant time losses due 
to flaring periods (particularly in PN) and in addition are now 
found to have strong residual quiescent contamination (> 30%). 
There is no cluster in the present sample found in these fields. 
The PN camera in field F32 is also significantly affected (~ 39% 
background enhancement). We identified two clusters (ID 476 
and 139) in this observation. Since the results from the PN and 
MOS cameras for these sources are in agreement within the er- 
ror bars we conclude that our background model was able to 
account for the enhancement. For more details on these sources 
see Sect. lA.2l 

The double component background model (see Sect. 13.1b 
used for the source detection and characterization can in prin- 
ciple account to first order for such an enhanced background 
by increasing the normalization of the background model. The 
vignetting function of such particle background has a different 
shape than the vignetting of the X-ray photons, but it is known 
only tentatively. We expect the errors from such first order ap- 
proximation to be small compared to other sources of uncer- 
tainty (including the shot noise). We thus decide to include into 
our analysis also fields with a strong residual quiescent contam- 
ination, but parameters derived for sources in these fields should 
be handled with caution. 

We treat out-of-time-events (OOTE) for the PN detector in 
a standard way. For each observation, we generate an OOTE 
eventlist with the epchain and remove time periods identi- 
fied in the two step cleaning process of the main PN eventlist. 
Whenever an image is extracted from the PN eventlist, we ex- 
tract also an image with the same selection criteria from the 
OOTE eventlist, scale this image with a factor of 0.063 (full 
frame readout mode) and subtract it from the PN image. 

3. Source detection 

As the main source detection algorithm we utilize the sliding 
box technique and a maximum likelihood source fitting in their 
improved implementation in the SAS tasks eboxdetect and 

- Slight contamination means < 15 % increase in the background with 
respect to normal levels (for details see lDe Luca & Molendi 2004) . This 
enhancement is modeled in first approximation by our composite back- 
ground model (see text further and Sect. 13. lb . 
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Fig. 1. Mosaic X-ray image of the 14 deg^ XMM-Newton sky survey. The false color image was constructed from the surface 
brightness in the 0.3 - 0.5, 0.5 - 2.0 and 2.0 - 4.5 keV bands. The white region (F) marks the 6 deg^ core of the survey presented 
in this work. Regions A, B and C constitute the extension of the survey by mosaic mode observations. The missing fields have 
significant losses due to soft proton flares. Bluer fields are affected by enhanced background. Green circles mark the positions of 
the present cluster sample and have a radius equal to rsoo. 



emldetect. A detailed description of the work flow and con- 
figuration of our detection pipeline, originally developed for the 
X MM-Newton Dista nt Cl uster Project (XDCP), c an be found 
in iFassbended (l2008h and iFassbender et al.l (l201 lal submitted); 
here we only summarize the main steps. 

Source detection is carried out in three different schemes: 

(i) the standard three band scheme: provides continuous, non- 
overlapping coverage in three energy bands: 0.3 - 0.5 keV, 
0.5 - 2.0 keV and 2.0 - 4.5 keV. 

(ii) the optimized single band scheme: covers the 0.35 - 2.4 keV 
band and was chosen to maximize the signal-to-noise-ratio for 
clusters of gal axies with a l arge range of redshifts and tempera- 
tures (see also ISchcjll2002h . This bandpass is expected to max- 



imize the signal-to-noise-ratio especially for high redshift sys- 
tems (z > 1) 

(iii) the five band spectral matched scheme: uses five partially 
overlapping bands (0.3 - 0.5, 0.35 - 2.4, 0.5 - 2.0, 2.0 - 4.5 and 
0.5 - 7.5 keV). This scheme is equivalent to a single band detec- 
tion in the full 0.3 - 7.5 keV range, where the energy intervals 
in the overlaps have higher weighting. The shape of the weight- 
ing function roughly m imics the expecte d continuum spectrum 
shape of a hot cluster (IFassbender! l2008h . This setup was used 
only to confirm detections from the first two schemes and we do 
not use any results derived from it in the current work. 
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Fig. 2. The 0.5-2 ke V band signal-to-noise ratio map of the XMM-BCS core region (region F in Fig. [T) smoothed with a Gaussian 
of 32" width. Circles indicate the rsoo radii of the detected clusters and are labeled with the cluster ID number in the catalog. Point 
sources have been subtracted using the method of .Finoguenov et aL (.20091) . 



3.1. Source list generation 

In order to obtain the raw source lists, we extract images from the 
cleaned eventlist for each detector and each band required in the 
given detection scheme (e.g. in the three band scheme three im- 
ages for each detector, in total nine images per field). We run the 
sliding box detection algorithm (eboxdetect in the so-called 
local mode) on these images. The background for each poten- 
tial source is estimated only locally in a detection cell of 5 x 5 
pixels in 4 successive runs with the number of pixels per cell 
doubled in each iteration. Sources detected by this procedure are 
then excised, creating an image usable for proper background 
estimation. 

We model the background of each detector and band individ- 
ually with a double component background model. This back- 
ground model is a linear combination of two templates based on 
vignetted and unvignetted exposure maps, taking into account 



the sky X-ray background (vignetted component) and the par- 
ticle and instrumental background (unvignetted in the first ap- 
proximation). 

The final sliding box detection is then run utilizing the fitted 
background model instead of a locally estimated background. 
For all sources above the detection threshold we carry out a max- 
imu m likelihood fitting (with the emld etect task). A beta pro- 
file (ICavaliere & Fusco-Femianolll976 l) with a fixed beta value 
of P — 2/3 convolved with the two dimensional point-spread 
function (PSF) is fitted to each source. The fit is carried out for 
all three detectors and all the bands in the given detection scheme 
simultaneously. The free parameters of the fit are the source po- 
sition, normalization of the model (for each detector and band) 
and the core radius, ^c, characterizing the source extent. If the 
extent of the source is not statistically significant, the source is 
refitted as a point source with extent fixed to zero. 
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Table 1. The individual XMM-Newton pointings. Quoted exposures are effective exposures with high background periods filtered 
out. 



Field ID RA (J2000) Dec (J2000) Exposure times (ks) 

OBSID Internal PN MOSl MOS2 





FOl 


23: 


:21 


:38 


4 


-56: 


:07: 


:34, 


4 


observation lost due to 


fl an Ti CT 

IIHL iiig 


1-/ J jT-^U 1 \J\J 1 


FOlb 


23: 


:22: 


00 


1 


-56: 


:09: 


:03. 


3 


7.8 


10.4 


10.4 




V/ J v/ J 1 


F02 


23: 


:24: 


:23 


.5 


-56: 


:07: 


13 


2 


observation lost due to 


fl ann c 

llaL lii^ 




F02b 


23: 


:24: 


:43 


.8 


-56: 


:09: 


03 





10.2 


13.2 


13.2 






F03 


23: 


:27 


:07 


.0 


-56: 


:07: 


15 


3 


observation lost due to 


fl arm tr 




F04 


23: 


:29 


:50, 


.6 


-56: 


:07: 


15 





5.2 


6.9 


6.9 






F04b 


23: 


30: 


: 1 1 


7 


-56: 


:09: 


01, 


2 


3.1 


12.7 


12.7 




VJ J V/J J O VJ J U 1 


F05 


23: 


:32 


:34, 


.6 


-56: 


:07: 


:12, 


.8 


observation lost due to 


fl ann c 




F06 


23: 


:35 


:39, 


.3 


-56: 


:08: 


18 


7 


5.6 


10.6 


10.6 






F07 


23: 


:20: 


:49 


.3 


-55 


:45: 


:35, 


1 


4.4 


9.6 


9.6 




(_/ J V/ J JOUOVJ 1 


F08 


23: 


:23 


:31 


4 


-55 


:45: 


:39. 


2 


9.3 


11.0 


9.8 






F09 


23: 


:26: 


:12, 


,7 


-55 


:46: 


10, 


,2 


2.3 


6.0 


6.0 






F09b 


23: 


26: 


: 1 1 


.6 


-55 


:46: 


:30. 


,2 


7.3 


9.8 


9.8 




v/»-' .JO J- \J\J X 


FIO 


23: 


:28: 


:55 


.3 


-55 


:45: 


:39, 


,2 


9.7 


12.6 


12.6 






Fl 1 


23: 


:31 


:37 


.8 


-55 


:45: 


:39. 


7 


7.2 


9.7 


9.7 






F12 


23: 


:34: 


:19, 


.5 


-55 


:45: 


:42, 


5 


10.8 


13.5 


13.5 






F13 


23: 


:37 


:01 


4 


-55 


:45: 


:39. 


,2 


2.3 


10.6 


10.6 




VI ^ .J O ± "VJ i- 


F14 


23: 


19 


:29, 


9 


-55 


:23: 


:01, 


1 


10.8 


13.9 


13.9 




O'SO'SSSI 'SOI 


F15 


23: 


:22 


:09 


7 


-55 


:23: 


:23, 


1 


7.4 


9.9 


9.9 






F16 


23: 


:24: 


:50, 


.3 


-55 


:23: 


:26, 


,3 


3.2 


11.7 


11.7 






F17 


23: 


:27 


:29, 


7 


-55 


:23: 


:45, 


9 


7.3 


10.0 


10.0 




OSO'S^SI 801 


F18 


23: 


:30: 
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The detection likelihood of a source is given by the detjnl 
parameter in the eboxdetect and emldetect tasks, defined as 
det_ml = -InPiand, where P,and is the probability of observed 
counts arising from pure random Poissonian fluctuations. In each 
step of the detection process, the minimum detection likelihood 
is set to 6, roughly equivalent to a > 3cr detection in terms of 
signal-to-noise ratio. 

The extent likelihood extjnl, defined analogously to char- 
acterize the probability of the source being extended, is required 
to be > 3 in the three-band scheme and > 5 in the single band 
scheme (corresponding approximately to minimum extent sig- 
nificances of ~ 2cr and ~ 3cr, respectively). 

For a more detailed discussion and justifi cation of the chosen 
detection schemes and thresholds we refer to lFassbendeil (l2008h . 



who also demonstrates the performance of the described source 
detection methods on over 450 archival XMM-Newton observa- 
tions in the framework of the XDCP project. A description of the 
used S AS tasks can be found in the S AS 7.1.0 reference manualQ 

In the current work, we aim for the best possible survey com- 
pleteness including the high redshift end of the cluster distri- 
bution and reliable source classification especially close to the 
detection thresholds. This is also helped by combining differ- 
ent detection schemes and setting relatively low extent thresh- 
olds. The increasing source contamination close to the detection 
threshold is treated with careful screening using the optical data 
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Fig. 3. Left: The black line shows the 10 second-binned PN light curve in the 0.3 - 10 keV band for the field F04. The beginning of 
the observation was affected by flaring. The green curve shows the light curve after the two-step cleaning (see Sect.|2]i, which safely 
removed all contaminated time intervals. Right: Examples of the detection pipeline products for field F04 in the 0.5 - 2 keV band 
of the PN detector: a) counts image, b) double-component background model, c) binary detection mask, d) reconstruction of all the 
detected sources. The green circle (2 arcmin radius) marks the cluster ID 018. 



and ancillary X-ray information (e.g. quality flags described in 
AppendixlAll. 

The detected sources create a raw master list of extended 
source candidates. Each of these candidates is then screened vi- 
sually with optical imaging data (4 band BCS imaging) and ac- 
cepted to the final cluster catalog only if a significant overdensity 
of galaxies in the photometric redshift space is found (Sect. |4]i. 
The available Spitzer imaging for the whole field will be used in 
the future to confirm z >1 systems, where the depth of the BCS 
imaging is not sufficient anymore. 

The purely X-ray based selection function will be devel- 
oped in subsequent work based on simulations, where complete- 
ness and contamination of different detection schemes will be 
studied. Guided by extensive simulations of X-ray observations 
jMiihlegger 2010.) . we will get a high precision description of 
the survey selection function. A statistically well defined cluster 
sample will be drawn from the current catalog (plus its 8 deg^ 
extension) and used to study the evolution of the cluster X-ray 
luminosity function and perform cosmological tests. 

3.1 .1 . Treatment of MOS CCDs in anomalous state 

A special note is required concerning the anomalous states of 
CCD#4 of the MOSl and CCD#5 of the MOS2 detectors and 
their effect on extended source detections. Half of our fields have 
the MOS2 CCD#5 in the anomalous state and ~ 20% have an 
anomalous MOSl CCD#4 (some observations are affected by 
both). These anomalous ("hot") states are characterized by high 
overall background count rates with atypical hardness ratios. The 
most affected are the soft bands (in our case 0.3 - 0.5 keV and 
0.5 - 2 keV bands, see also Kuntz & Snowden 2008). 

We check for the presence of a hot chip in an observa- 
tion by comparing count rates extracted from the suspected 



chip and the mean count rate of three other chips in symmet- 
rical positions around the central chip (i.e. the mean count rate 
of CCD#2, CCD#4, CCD#7 of MOS2 and CCD#3, CCD#5, 
CCD#7 of MOSl). These reference chips were selected, because 
they best match the area, shape and position of the affected chips 
(see middle panel of Fig. |4|i. The count rates calculated in the 
0.3 - 2.0 keV band from the three reference chips are then aver- 
aged to reduce shot noise and a chip is flagged hot, if its count 
rate is more than 10% higher than the mean count rate from the 
reference chips. This threshold is chosen to be very conservative 
and was found to perform excellently, since chips in anomalous 
states have typically count rates 50 - 100% higher than the ref- 
erence rate. The algorithm also automatically flags observations 
where a bright source lies on a reference chip. In these cases we 
manually excise the source and repeat the calculation to obtain 
an unaffected count rate ratio. 

The exceptionally high background of the hot chips leads to 
many spurious extended source detections, when left untreated 
(see Fig. |4]i. We flag sources as possibly spurious detections 
caused by the presence of a hot chip if at the same time: 1) they 
lie on a chip that was flagged hot, 2) are extended, 3) the de- 
tection likelihood from the given hot MOS detector in the soft 
bands (sum of the 0.3 - 0.5 and 0.5 - 2.0 keV bands) accounts 
for more than 90% of the total detection likelihood and 4) the 
source would be under our detection threshold without the detec- 
tion on the affected chip. We still visually checked every ffagged 
source also in the optical images and confirmed the classification 
of these sources as spurious. 

An example of this procedure can be seen in Fig.|4] The ob- 
servation of field F21 has a hot MOS2 CCD#5, clearly visible 
as an enhanced background (in the raw image in the left panel 
and in the model background in the middle). The 8 extended 
sources detected on this chip were flagged as potentially spu- 
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Fig. 4. Left: Image of field F2 1 taken by the M0S2 camera in the 0.5-2.0 keV band. The MOS2 CCD#5 is visibly in an anomalously 
high ("hot") state with an enhanced background. Sources detected in this field are marked by green circles. Sources with red circles 
were automatically flagged as possibly spurious detections caused by the presence of the hot chip. Middle: A composite background 
model for the same detector and band created by fitting the double component model independently to the CCD#5 and the rest of 
the chips. The three blue-marked chips are the reference chips used to identify hot chips in the observations. Right: The ratio of 
the total detection likelihood (log scale) from the MOS2 CCD#5 in the 0.3 - 0.5 and 0.5 - 2.0 keV bands to the total detection 
likelihoods from all other detectors and bands (log scale). Blue bars show the confirmed clusters from our sample, the red bars 
the 8 flagged sources from field F21 (from the left panel). The vertical line marks where the soft band MOS2 detection constitutes 
90% of the total detection likelihoods in all detectors. The flagged sources were confirmed as spurious by the optical data. A single 
confirmed cluster (ID 275) appears above the threshold, but is not flagged as spurious since it would have been above the detection 
likelihood even without the MOS2 detection (i.e. not meeting all the required criteria described in Sect. l3.1.TT i. 



rious based on the described criteria. The detection likelihood 
ratio (the MOS 2 detection likelihood in the soft bands over the 
total detection likelihood) of these 8 sources are displayed on the 
left panel of Fig. |4](red) as compared to the sample of confirmed 
clusters in our sample (blue). 

A similar criterion can be applied in principle also to spu- 
rious point source detections. An additional improvement can 
be achieved by weighting the input detection likelihoods by the 
number of pixels in the detection aperture in order to avoid a 
possible bias, if a source has a low detection likelihood in one of 
the reference detectors only because it falls on a chip gap or is 
(partially) out of the field-of-view. 

We also make an attempt to model the high background of 
the hot chips by fitting the double component model to a hot chip 
(in first approximation) and another double component model to 
the remaining chips. The two parts of the background model are 
then combined to create a composite background map for the full 
detector area (middle panel of Fig. |4]i. All the extended sources 
on hot chips flagged as spurious with the described detection 
likelihood test, are not detected when the composite background 
maps are utilized, confirming the reliability of our classification. 
The effect of using a composite background instead of a stan- 
dard background on detections coming from the remaining, non- 
anomalous chips is minor, since the two background models in 
these areas diff'er typically by less than 5%, and only the soft- 
est bands of each detection scheme are affected. For the source 
characterization in observations affected by hot chips we use ex- 
clusively composite background maps. 

3.2. Growth curve analysis 

The X-ray count rate is the most direct cluster observable. With 
an estimate of the energy conversion factor at hand (see also 
Sect. I3.3l l we can calculate from it the X-ray flux Fx, which in 
turn can be converted to X-ray luminosity Lx- The luminosity is 
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Fig. 5. Example of the growth curve analysis of source ID 018 
(photo-z=0.39). The cluster's redshift and luminosity are close 
to the median values of the entire sample. The curves show the 
encircled cumulative flux as a function of radius (PN: blue curve, 
combined MOS: red, dot-dashed). The PN and MOS curves are 
in good agreement. Dashed lines mark the flux measurement er- 
ror bars which include the Poisson noise and an additional 5% 
systematic error from the background estimation. The estimated 
plateau flux is F;,/„, = 5.19x10"'^ erg s"' cm"^ (horizontal line), 
reached at rpiat ~ 90 arcsec. The vertical line signifies the esti- 
mated rsoo radius of the source, rsoo = 0.6 Mpc (~ 117 arcsec). 
In this case, the plateau radius is slightly smaller than rsoo and 
the flux and luminosity for rjoo had to be extrapolated from their 
plateau values. The required extrapolation is only ~ 2% in this 
case. See Sect. l3.3l for details. 
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a key parameters since it allows us to calculate other important 
physical parameters (particularly the mass of the system) from 
scaling relations. 

A typical cluster of galaxies in relatively shallow observa- 
tions like ours appears as a faint diffuse source with typically 
> 100 source photons registered (total from all three standard 
bands and detectors). Thus in order to get a reliable measure- 
ment of the flux and trace the emission of the cluster as far out 
as possible, we have to employ a robust method. In this work we 
utilize the growth curve method developed for the REFLEX and 
N ORAS cluster samples d erived from the ROSAT all sky survey 
bv lBohringer et al] (1200 0). Here we summarize the procedure. 

For each source, we extract images, exposure maps and 
background maps in the 0.5 - 2.0 keV band, excluding all point 
sources detected by the pipeline. MOS 1 and MOS2 products are 
then directly co-added, since the difference in their response ma- 
trices is small. We run the growth curve program on the PN and 
co-added MOS images independently. 

In this analysis we use the X-ray center coordinates obtained 
from the beta model fitting procedure in the source detection 
step. We also explored the possibility of recentering by min- 
imizing the dipole m oment of the count distribution (see e.g. 
iBohringer et al.ll20I0l) . This procedure usually yielded centers 
very close to the best-fit beta model coordinates, but for faint 
sources often completely diverged. The best-fit coordinates were 
always found to be a good description of the detected X-ray 
emission centroid. 

Counts are extracted from the image in concentric rings start- 
ing from the center and scaled by the exposure time. In this 
way we obtain the total (source + background) count rate pro- 
file. The expected background count rate is estimated from the 
model background map and subtracted for each ring from the to- 
tal count rate, obtaining the source count rate profile. The growth 
curve is the cumulative background subtracted source count rate 
profile. An example of a growth curve is displayed in Fig. |5] 
(shown here in flux units using the energy conversion factor cal- 
culated as described in Sect. 13.3b . 

We term the radius of the full aperture inside which a stable 
growth curve can be obtained, the extraction radius rgxt (typi- 
cally 150 - 200 arcsec). It is adjusted for each source individu- 
ally (increased for brightest, most extended sources or trimmed 
for sources close to the edge of FOV or to a partially blended 
systems) and includes the source itself as well as enough sky 
region to check the reliability of the double component back- 
ground subtraction. 

If the background model describes the local background ac- 
curately, the growth curve levels off" to a flat plateau at the outer 
edge of the source. To estimate the total detected cluster emis- 
sion, we first calculate the significance radius rsig, defined as the 
radius outside which the source signal increases less than the Icr 
uncertainty in the count rate. The significance radius thus gives 
the outermost radius where the potential increase of the growth 
curve becomes less than Icr significant. 

To alleviate the effect of shot noise, rsig is determined by 
smoothing the growth curve in 20 and 48 arcsec windows (5 and 
12 pixels respectively). For most clusters the two estimates are in 
agreement. In the remaining cases, the local background usually 
exhibits irregular features not captured by the double component 
model and we select the more appropriate rsig and plateau after 
visual inspection. 

In addition, a single multiplicative correction factor to the 
background model can be set, if the plateau exhibits a signif- 
icant residual slope. This additional factor corrects the over- 
all normalization of the double component model locally in- 



side rext- The average background correction factors are -2% 
(i.e. a 2% decrease compared to the default double component 
background) for PN and 0% for MOS (with standard deviations 
7% and 8%, respectively). More than 3/4 of t he present sam- 
ple have correction factor s smal ler than 10%. IReiprichI (1200 Ih 
and lReiprich & Bohringerl (|2002|) used a similar correction pro- 
cedure utilizing a second order polynomial to obtain stable 
plateaus. In our case, a simple correction factor turned out to 
be sufficient and not leading to background over-fitting. 

After setting the background correction, the total source 
count rate is estimated as the count rate of the plateau. The flat 
plateau of the growth curve outside rsig is then fitted with a line. 
If the slope of the line is less than 0.8% per radial bin, the plateau 
fit is accepted and the plateau count rate CTRpiat is estimated as 
the mean of the fitted line. If the slope is still not negligible, 
an additional attempt is made to find a stable plateau by itera- 
tively removing the outermost and innermost (still outside rsig) 
bins. We note that in ~ 80% of cases the first simple fit is fully 
acceptable and no further iterations are necessary. For more de- 
tailed description of the iterative process and quality flags of the 
plateau fit see Sect. |A]in the appendix. The plateau radius rpiat, 
is defined simply as the radial distance where the growth curve 
first reaches CTRpiat- 

We provide a performance test of our X-ray photome- 
try method on the example of the XMM-LSS cluster catalog 
(iPacaud et al.ll2007 l) in Appendix IC. 21 The main advantages of 
the growth curve method thus are: (i) Excellent sensitivity al- 
lowing us to trace cluster emission to the outermost faint out- 
skirts, (ii) It makes no assumptions about the source profile un- 
like methods based on beta model fitting, which is fully degener- 
ate in the regime with < 400 - 500 counts and is known not to be 
an appropriate description of cluster emission for irregular and 
cool core clusters, (iii) The method allows to check and correct 
the background modelling which is done for the whole field of 
view, by adjusting several parameters to the conditions local to 
each analyzed source, (iv) The PN and combined MOS growth 
curves are treated completely independently. Their comparison 
provides us with an important consistency check and allows us to 
treat instrument specific features in the background separately. 

3.3. Physical parameter estimation 

With a stable PN and MOS growth curve at hand we determine 
all the relevant physical parameters of the clusters (see Table |2]) 
in an iterative way: 

1) The physical parameters are set to their initial values 
(r50o,r5oo) = (0.5Mpc,2.5keV). 

2) The physical aperture radius is converted to arcseconds 
using the assumed cosmological parameters. The total count rate 
inside this radius is estimated from the PN and MOS growth 
curves in the 0.5 - 2 keV band. 

3) The count rates are converted to flux with an energy 
conversion factor (ECF) calculated assuming a MEKAL spec- 
ti-al model (iMewe et al.lll985tlKaastrallI992l:lLiedahl et al.lll995h 
with abundance set to 0.3 times solar abundance, temperature 
equal to the trial r^oo value and redshift set to the photo-z value 
(or spectroscopic redshift where available). 

To account for the spatial variation of the spectral response 
of the detectors we calculate a response matrix for each source 
individually in a 150 arcsec aperture centered on the source for 
the THIN filter The MOS2 response matrix is used to calculate 
the ECF of the co-added MOS count rates. 
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This procedure is used to estimate the 0.5 - 2 keV and bolo- 
metric fluxes for PN and MOS and the same model is also used 
to convert the fluxes to luminosities. 

4) In some cases the estimated value of rgoo in the given iter- 
ation is larger than rpiat (i.e. larger than the region with directly 
measurable emission), and therefore an extrapolation factor has 
to be applied to the flux and luminosity estimates. We correct for 
the missing flux between the plateau radius and current iteration 
estimate of rgoo by extrapolating the source emission assuming a 
beta model. The /? and rcore parameters of the beta model are cal- 
culated using the scaUng re l ation s of Reiprich (2001) (see also 
iReiprich & Bohringedl2002l; iFinoguenov et al.ll2007h : 




(1) 
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Table 2. Physical parameters of the clusters sample. Fluxes and luminosities are given for the 0.5 - 2 keV band. See Appendix lA.2l for notes on individual sources. Notes: ' the cluster was detected 
in observations strongly affected by flaring; " the cluster is heavily affected by blending with a nearby source. Notes for redshifts: ° spectroscopic redshift; '' a high redshift system for which a secure 
photometric redshift estimate is not possible from the current photometric catalog (the listed parameters are tentative estimates). 



ID 


R.A. (J2000) 


Dec (J2000) 


z 


/"pl;U 


^ plal 


csoo 


F50Q 


^-500 


7^500 


M500 


1^500 


M200 




(deg) 


(deg) 


(photo.) 


(ai-cmin/r^fli,) 


(10"" erg s"' cm"-) 


(kpc) 


(10 '* erg s 'cm -) 


(10"" erg s"') 


(keV) 


(10'-' Mo) 


(10" Mo keV) 


(10" Mo) 


011^ 


351.8070 


-56.0615 


0.97 ±0.10 


0.7/0.6 


2.80 ± 0.42 


577 ± 54 


2.85 ± 0.43 


11.8 ± 1.8 


3.4 ± 1.0 


16.4 ±4.6 


4.8 ±2.9 


24.3 ± 6.8 


018 


352.4828 


-56.1360 


0.39 ± 0.04 


1.5/0.7 


5. 14 ±0.50 


633 ±58 


5.21 ±0.51 


2.6 ± 0.3 


2.3 ±0.7 


10.9 ±3.0 


1.8± 1.1 


15.2 ±4.2 


032 


352.1778 


-55.5662 


0.83 ± 0.07 


1.6/1.0 


8.06 ± 0.68 


702 ± 64 


8.00 ± 0.66 


21.5 ± 1.8 


4.2 ± 1.2 


25.0 ±6.9 


10.8 ±6.4 


37.0 ± 10.2 


033 


352.0448 


-55.8400 


0.79 ± 0.05 


1.1/0.9 


2.94 ± 0.25 


593 ± 54 


2.94 ± 0.25 


7.6 ± 0.6 


3.0 ±0.9 


14.4 ± 4.0 


3.5 ±2.1 


2I.0±5.8 


034 


352.6538 


-55.7270 


0.28 ± 0.02 


1.7/0.8 


2.84 ± 0.45 


536 ±50 


2.89 ± 0.46 


0.7 ±0.1 


1.5 ±0.4 


5.8 ± 1.6 


0.5 ±0.3 


7.9 ±2.2 


035 


353.4388 


-55.6387 


0.67 ± 0.05 


1.2/0.9 


2.02 ± 0.32 


560 ± 53 


1.95 ±0.31 


3.6 ± 0.6 


2.4 ±0.7 


10.5 ±3.0 


1.8± 1.1 


15.0 ±4.2 


038" 


353.5130 


-55.8156 


0.39 ± 0.05 


1.1/0.7 


1.30 ±0.20 


503 ± 47 


1.33 ±0.20 


0.7 ±0.1 


1.5 ±0.4 


5.4 ± 1.5 


0.4 ±0.3 


7.5 ±2.1 


039 


349.8214 


-55.3244 


0.18 ±0.04 


2.3/0.7 


9.12 ±0.56 


586 ±53 


9.33 ± 0.57 


0.8 ±0.1 


1.6 ±0.5 


6.8 ± 1.9 


0.7 ±0.4 


9.2 ±2.5 


044 


349.2212 


-54.9036 


0.44 ± 0.02 


2.4/1.1 


17.14± 1.15 


787 ± 72 


16.75 ±1.10 


10.5 ±0.7 


3.6 ± 1.0 


22.2 ±6.1 


7.5 ±4.4 


31.6 ±8.6 


069 


350.9631 


-54.8923 


0.75 ± 0.07 


1.7/1.3 


2.22 ± 0.43 


560 ± 53 


2.03 ± 0.35 


4.8 ±0.8 


2.6 ± 0.7 


1I.6±3.3 


2.2 ± 1.4 


16.7 ±4.7 


070 


350.6286 


-54.2691 


0.152° 


2.9/0.6 


80.78 ± 1.64 


815 ±74 


83.43 ± 1.70 


5.0 ±0.1 


2.9 ±0.8 


17.9 ±4.9 


4.5 ±2.6 


24.6 ± 6.7 


081 


351.8470 


-55.2624 


0.85 ±0.12 


0.4/0.4 


1.01 ±0.17 


497 ± 47 


1.05 ±0.18 


3.5 ±0.6 


2.3 ±0.7 


9.1 ±2.6 


1.4 ±0.9 


13.2 ±3.7 


082 


351.5779 


-55.3859 


0.63 ± 0.05 


0.9/0.7 


1.32 ±0.23 


526 ± 50 


1.34 ±0.23 


2.2 ± 0.4 


2.0 ±0.6 


8.3 ±2.4 


1.1 ±0.7 


11.8 ±3.3 


088 


352.1748 


-55.2234 


0.43 ± 0.04 


2.5/1.3 


8.50 ± 1.27 


679 ± 63 


7.43 ± 1.01 


4.6 ± 0.6 


2.7 ±0.8 


14.1 ± 3.9 


3.0 ± 1.8 


19.8 ± 5.5 


090 


352.2366 


-55.4081 


0.58 ± 0.02 


0.8/0.7 


0.76 ±0.18 


481 ±47 


0.78 ±0.18 


1.1 ±0.2 


1.6 ±0.5 


6.0 ± 1.8 


0.6 ± 0.4 


8.3 ±2.4 


094 


353.0185 


-55.2120 


0.269" 


1.3/0.6 


4.83 ± 0.53 


583 ± 54 


5.02 ± 0.55 


1.1 ± 0.1 


1.7 ±0.5 


7.4 ±2.1 


0.8 ±0.5 


10.1 ± 2.8 


109 


351.9058 


-54.2705 


1.02* 


0.7/0.7 


1.45 ±0.21 


510 ±48 


1.47 ±0.21 


7.3 ± 1.1 


2.8 ±0.8 


12.0 ±3.4 


2.6 ± 1.6 


17.7 ± 5.0 


110 


352.5161 


-54.2388 


0.47 ± 0.06 


1.4/0.8 


3.77 ± 0.49 


605 ± 56 


3.47 ± 0.47 


2.7 ± 0.4 


2.3 ±0.7 


10.4 ±2.9 


1.7 ± 1.0 


14.6 ±4.1 


126 


351.6393 


-55.0206 


0.42 ± 0.02 


1.9/1.0 


5.48 ±0.41 


643 ± 59 


5.38 ± 0.40 


3.2 ± 0.2 


2.4 ± 0.7 


11.8 ±3.2 


2.1 ± 1.3 


16.5 ±4.5 


127 


351.8492 


-55.0648 


0.207° 


0.8/0.3 


2.58 ±0.29 


499 ± 46 


2.91 ±0.33 


0.4 ±0.1 


1.2 ±0.4 


4.4 ± 1.2 


0.3 ± 0.2 


5.9+ 1.6 


132 


352.0084 


-54.9292 


0.96 ±0.17 


1.3/1.1 


2.91 ±0.34 


571 ±53 


2.74 ±0.31 


11.1 ± 1.2 


3.3 ±0.9 


15.7 ±4.4 


4.4 ±2.7 


23.3 ± 6.5 


136 


350.5036 


-54.7500 


0.36 ± 0.02 


2.4/1.1 


8.29 ± 0.76 


673 + 62 


8.15 ±0.68 


3.4 ±0.3 


2.5 ± 0.7 


12.6 ±3.5 


2.4 ± 1.4 


17.6 ±4.8 


139+° 


351.3953 


-54.7212 


0.169° 


1.9/0.6 


4.83 ± 0.65 


517±48 


5.01 ± 0.67 


0.4 ±0.1 


1.3 ±0.4 


4.6 ± 1.3 


0.3 ±0.2 


6.2 ± 1.7 


150 


352.5015 


-54.6184 


0.176° 


2.9/0.8 


18.25 ±0.89 


654 ± 59 


18.48 ±0.91 


1.6 ±0.1 


2.0 ± 0.6 


9.4 ±2.6 


1.3 ±0.7 


12.9 ±3.5 


152 


352.4168 


-54.7886 


0.139° 


0.9/0.3 


2.50 ±0.41 


448 ±42 


2.97 ± 0.49 


0.2 ±0.1 


1.0 ±0.3 


2.9 ±0.8 


0.1 ±0.1 


3.9± 1.1 


156 


353.8815 


-54.5865 


0.67 ± 0.06 


0.8/0.6 


3.35 ± 0.24 


614 ±56 


3.43 ± 0.24 


6.0 ±0.4 


2.8 ± 0.8 


13.9 ±3.8 


3.2 ± 1.9 


19.9 ±5.5 


158 


353.6032 


-54.4586 


0.55 ± 0.03 


1.4/0.9 


3.57 ±0.54 


617 ±58 


3.58 ±0.54 


4.0 ±0.6 


2.5 ± 0.7 


12.1 ±3.4 


2.3 ± 1.4 


17.2 ±4.8 


210 


353.5240 


-55.7859 


0.83 ± 0.09 


0.7/0.7 


0.55 ±0.1 1 


451 ±43 


0.56 ±0.11 


1.9 ±0.4 


1.9 ±0.5 


6.6 ± 1.9 


0.8 ± 0.5 


9.5 ±2.7 


227 


350.5425 


-55.4199 


0.346° 


1.1/0.6 


1.51 ±0.19 


506 ± 47 


1.56 ±0.20 


0.6 ±0.1 


1.4 ±0.4 


5.3 ±1.5 


0.4 ±0.3 


7.2 ±2.0 


245 


351.0160 


-55.0225 


0.62 ± 0.03 


0.8/0.7 


0.97 ± 0.18 


500 ± 48 


0.99 ± 0.18 


1.6 ±0.3 


1.8 ±0.5 


7.0 ±2.0 


0.8 ± 0.5 


9.9 ±2.8 


275 


353.6991 


-55.2736 


0.29 ± 0.03 


1.7/0.8 


2.93 ± 0.45 


541 ±51 


2.95 ± 0.45 


0.8 ±0.1 


1.5 ±0.4 


6.0 ± 1.7 


0.5 ±0.3 


8.2 ±2.3 


287 


354.2119 


-55.2988 


0.57 ± 0.04 


0.9/0.7 


0.98 ±0.35 


501 ±54 


1.00 ±0.36 


1.3 ±0.5 


1.7 ±0.5 


6.7 ±2.1 


0.7 ±0.5 


9.3 ±3.0 


288 


353.7523 


-54.9164 


0.60 ± 0.04 


1.8/1.2 


2.75 ± 0.63 


580 ±56 


2.44 ± 0.48 


3.4 ±0.7 


2.4 ± 0.7 


10.7 ±3.1 


1.9 ± 1.2 


15.3 ±4.4 


357 


353.6200 


-54.6066 


0.48 ± 0.06 


2.0/1.2 


3.56 ±0.47 


596 ± 55 


3.16 ±0.40 


2.6 ±0.3 


2.2 ±0.6 


10.1 ±2.8 


1.6 ±0.9 


14.2 ±3.9 


386 


353.9763 


-56.0928 


0.53 ± 0.05 


0.7/0.6 


0.66 ±0.17 


468 ± 47 


0.67 ±0.18 


0.8 ± 0.2 


1.5 ±0.4 


5.2 ± 1.6 


0.4 ±0.3 


7.2 ±2.2 


430 


351.3891 


-55.7327 


0.206° 


1.3/0.6 


1.59 ±0.33 


455 ± 44 


1.69 ±0.36 


0.2 ±0.1 


1.0 ±0.3 


3.3 ± 1.0 


0.2 + 0.1 


4.4 ± 1.3 


444 


354.0839 


-55.5189 


0.71 ±0.05 


0.9/0.7 


1.27 ±0.25 


521 ±50 


1.29 ±0.26 


2.8 ± 0.6 


2.2 ±0.6 


8.8 ±2.5 


1.3 ±0.8 


12.6 ±3.6 


457 


352.1177 


-54.2472 


0.1° 


1.6/0.5 


1.98 ±0.52 


387 ±39 


2.29 ± 0.60 


0.1 ±0.1 


0.7 ±0.2 


1.8 ±0.5 


< 0.1 


2.4 ± 0.7 


476+° 


351.4166 


-54.7412 


0.102° 


2.3/0.5 


10.27 ± 1.18 


508 ± 47 


11.04± 1.27 


0.3 ±0.1 


1.2 ±0.3 


4.1 ± 1.1 


0.2 ±0.1 


5.5 ±1.5 


502 


349.9334 


-54.6400 


0.55 ± 0.05 


0.7/0.5 


1.55 ±0.14 


540 ± 50 


1.61 ±0.14 


1.9 ±0.2 


2.0 ± 0.6 


8.1 ±2.2 


1.1 ±0.6 


11.4 ±3.2 


511 


353.0628 


-54.7006 


0.269° 


1.5/0.6 


4.41 ±0.71 


574 ± 54 


4.53 ± 0.73 


1.0 ±0.2 


1.7 ±0.5 


7.1 ±2.0 


0.7 ± 0.4 


9.7 ±2.7 


527 


350.5456 


-56.3127 


0.79 ± 0.06 


1.0/0.8 


2.26 ±0.37 


568 ± 54 


2.27 ± 0.37 


6.0 ± 1.0 


2.7 ±0.8 


12.6 ±3.6 


2.7 ± 1.7 


18.3 ±5.2 


528 


353.8357 


-55.5442 


0.35 ± 0.02 


0.8/0.5 


0.63 ±0.21 


441 ±46 


0.68 ± 0.22 


0.3 ±0.1 


1.1 ±0.3 


3.5 ±1.1 


0.2 ±0.1 


4.8 ± 1.5 


538 


353.5258 


-54.7310 


0.20 ± 0.02 


1.6/0.7 


1.99 ±0.81 


469 ± 52 


2.07 ± 0.85 


0.2 ±0.1 


1.1 ±0.3 


3.6 ± 1.2 


0.2 + 0.1 


4.8 ± 1.6 


543 


353.1806 


-54.8297 


0.57 ± 0.03 


1.7/1.0 


4.15 ± 1.13 


629 ±61 


4.07 ± 0.90 


4.9 ± 1.1 


2.7 ±0.8 


13.2 ±3.8 


2.8 ± 1.7 


18.7 ±5.5 


547 


351.0815 


-55.4305 


0.241° 


1.1/0.5 


1.05 ±0.33 


443 ±46 


I.I2±0.35 


0.2 ±0.1 


1.0 ±0.3 


3.1 ± 1.0 


0. 1 ± 0. 1 


4.2 ± 1.3 
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The good sensitivity of the growth curve method allows us 
to trace cluster emission out to large radii, therefore the re- 
quired extrapolation is typically minor. The mean correction is 
~ 5% I ~ 6% for PN/MOS (~ 46% at maximum). 

In cases when rsoo < rpiat no extrapolation is needed and the 
flux and luminosity estimates are independent of any assumption 
about the spatial distribution of the source emission. 

5) The source flux and luminosity are then obtained by aver- 
aging the PN and MOS fluxes (luminosities) weighted by their 
inverse squared errors. Sources for which the PN and MOS esti- 
mates do not agree or one of the estimates is missing (e.g. source 
outside of the FOV of a given detector) are flagged (Table lA.ll 
in the appendix). An X-ray photometric quality flag is also as- 
signed to each source based on the quality of the plateau fit, por- 
tion of pixels outside the detection mask, presence of anomalous 
features in the X-ray background and visual screening. 

6) We then use this total (i.e. camera averaged) bolometric 
luminosity value to calculate the temperature and mass of the 
cluster for the nex t iteration by utili zing the L - T and L - M 



cluster tor tne nex t iteration by utui zi 
scaling relations of iPratt et al] ( 120091) : 



1.38 X 1044 erg s- 



M = 2xlO"*M(; 



1/2.08 



r = 5 keV 



7.13X 1044 erg s- 



1/3.35 



(2) 



(3) 



These relations were ob tained by the BCES orthogonal fits 
(lAkritas & Bershadvl[T996l) . which do not treat T^oo (or M500) 
as the independent variable, since our independent variable is in 
fact Lx(< rsoo). At this stage it is impossible to safely detect 
and remove emission from possible cooling cores because of the 
limited resolution of XMM-Newton. Therefore we opt not to do 
so and use the relations that include also the core regions. 

The F500 parameter is calculated from the .Pratt et al.i (l2009t) 
relation: 



Fx = 2x lO'^^MokeV 



5.35 X 1044 erg s- 



1/1.04 



(4) 



7) The M500 estimate is then used to obtain a new rjoo radius 
from r5oo = -^3M5oo/47r500pc(-z)^ where pciz) is the critical den- 
sity of the Universe at redshift z- This new rgoo aperture is used 
along the updated Tsqo value to recalculate the fluxes and lumi- 
nosities by repeating steps 1-7. The process is repeated until 
convergence to a final solution. 

The final parameters are listed in Table |2l The table also 
includes mass estim ates in the ^200 aperture. An NFW profile 
(Navarro et alj |1997h was assumed in order to extrapolate the 
mass from the rsoo to r2oo aperture. The extrapol ation factor was 
calcul ated using the approximation derived by iHu & Kravtsovl 
(l2003h . where the parameters of the NFW profile are iteratively 
est imated from the final (i.e. converged) M500 mass by utilizing 
the Bull ock et alj (1200 ll) relation for the concentration parameter 
calculation. 

We provide a discussion on the scaling relations utilized here 
and the error budget of our iterative method in Sect j6.1l 

4. Photometric redshift estimation 

In order to measure the photometric redshifts (photo-zs) of 
the X-ray selected systems in our sample, we applied the red- 
sequence redshift estimator to the Blanco Cosmology Survey 



(BCS) imaging data which covers two 50 deg^ patches of the 
southern sky and includes the full area of the present XMM- 
Newton survey. The BCS is a 60 nights NOAO survey program 
carried out from 2005-2008 on the Blanco 4m telescope at Cerro 
Tololo Interamerican Observatory. This griz survey was tuned to 
the required depths to follow the passively evolving cluster red 
sequence population to L* + 1 at 5cr significance out to z = 1. 
The data were acquired using two offset layers of imaging in g 
and r band, and three offset layers of imaging in / and z band. 

The BCS data have been processed and calibrated using a de- 
velopment version of t he Dark Energy Su rvey data management 
system (DESDM v4, iMohr et all |2008|) . The core processing 
includes flat fielding, bias subtraction, illumination and fringe 

corrections. Astrometric refinement uses the USNO-B2 catalog 

1 III I ^ 

and the AstrOmatic tool SCAMP (Bertin 2006). Reduced single 
epoch images are combined into deep, coadded images using the 
AstrOmatic tool S Warp together with DESDM code that homog- 
enized the PSF to the median seeing within each tile/band com- 
bination. Model fitting photometry using bulge-i-disk decompo- 
sition was carried out using the extended version of SExtractor 
(iBertin & Arnouts 199 6) developed within DESDM to enable 
PSF corrected model fitting to all detected objects on an image. 
Photometric calibration was carried out using PSF model magni- 
tudes calibrated using the stellar locus in the col or space defined 



by grizH where the H band data from 2MASS (Skrutskie et al 
2000) provides the overall photometric scale (.Armstrong et al 
1201 0). A similar approach has been used for the processing 
and calibration of all SPT cluster followup and redshift esti- 
mation using Blanco 4m data, and the results have enabl ed de- 
tailed studies of cluster galaxy properties (IZenteno et alj|201 Ih 
as well as precise photometr ic redshift estimation ( High et al] 
l20Tol:IWlliamson et al.ll20Tll) . 

The ful l description of ou r photometric redshift method is 
provided in lSong et al.l (1201 ll) : here we give its brief summary. 
The red-sequence redshift estimator utilizes all available filters, 
(g-, r-, i-, and z) to search for redshift peaks in the density dis- 
fiibution of galaxies within a radius of 0.8 Mpc centered on 
the X-ray detection. To define the red-sequence at each red- 
shi ft slice, we assume a sing le stellar population (SSP) model 
bv lBruzual & CharlotI (120031) with a single burst of star forma- 
tion at z/ = 3 and passive evolution of red galaxies thereafter. 
The SSP models are run with six distinct metallicities in order 
to be able to model a tilted red sequence. The models are cal- 
ibrated to reproduce t he tilt of the color -magnitude relation for 
the Coma cluster (e.g. lBower et aljri992h . 

A single stellar population (SSP) model assuming a single 
burst of star formatio n at Zf = 3 and passive e volution of red 
galaxies thereafter, by iBruzual & Chariot ( l2003h is used to de- 
fine the red-sequence at each redshift slice, which is calibrated 
to Coma cluster 

The contribution of background galaxies is estimated from a 
surrounding 36' x 36' sky patch and statistically subtracted. For 
each X-ray cluster candidate the whole redshift range from z - 
to z = 1 .05 is scanned through using simultaneously two colors 
that bracket the 4000 A break for a given redshift. This sup- 
presses false overdensity peaks at transitional redshifts where 
the 4000 A break moves between two adjacent bands (e.g. the 
transition between the g and r band around z ~ 0.35). Once a 
peak in redshift space is identified, we refine the redshift esti- 
mate by fitting a Gaussian function to the redshift density distri- 
bution. We then select cluster members in a stripe (0.05 width in 
color) around the estimated red-sequence. The final cluster red- 
shift value is calculated as the inverse color error weighted mean 
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Table 3. Spectroscopic redshift for 12 clusters in the redshift 
range z - - 0.4. The redshifts were obtained from long-slit 
spectroscopic observations at the NTT telescope. The redshifts 
of the BCG galaxies are in the z^^^ column. Four clusters have 
a redshift for one additional member galaxy (zfpec)- Photometric 
redshifts Zphoto are taken from Table |2l For five systems we also 
provide the photometric redshifts from the SCS survey (M09, 
MIO). 



ID 


^spec 


z" 

^spec 


■Spholo 


,scs 

^pholo 


70 


0.152 


0.152 


0.17 ±0.03 


0.12 


94 


0.269 




0.29 ± 0.04 




127 


0.207 


0.209 


0.22 ± 0.02 




139 


0.169 




0.18 ±0.01 




150 


0.176 


0.173 


0.20 ± 0.02 


0.14 


152 


0.139 




0.17 ±0.02 




227 


0.346 




0.35 ± 0.04 




430 


0.206 


0.205 


0.18 ±0.01 




457 


0.1 




0.10 ±0.01 




476 


0.102 




0.10 ±0.01 


0.1 


511 


0.269 




0.26 ± 0.02 


0.2 


547 


0.241 




0.22 ± 0.02 


0.18 



redshift of the selected member galaxies. This assures that the 
reliability of the photo-z values for the whole system is always 
better than for any individual galaxy. In a few cases two or more 
solutions were found by our algorithm. For these systems we 
visually check the obtained redshift distributions and select the 
more likely solution given the positions of galaxies with respect 
to the X-ray emission. An example of a final color-magnitude di- 
agram is shown in Fig.|6]for cluster ID 018 (its redshift is close 
to the median redshift of the cluster sample). 

The described photo-z estimation method allows us to mea- 
sure the cluster redshift with good precision up to z ^ 0.8 even 
for low richness systems. The overall uncertainty of the photo-zs 
is on the ~ 10% level. While care was taken to obtain reliable 
results also for z > 0.8 systems (see Fig. Qfor two examples), 
here the already obtained Spitzer mid-infrared observations will 
provide an important improvement in subsequent work. The fi- 
nal photometric redshifts are presented in Table |2] A more de- 
tailed analysis of optical counterparts for our systems including 
optical luminosity and richness estimates will be presented in a 
companion paper (Song et al., in prep.). 

4.1. Spectroscopic redshifts 

Spectroscopic redshifts are required to identify the clusters as 
compact objects, to derive precise physical parameters and later 
for cosmological modeling. In order to make a first step towards 
these goals we have carried out spectroscopic observations of a 
subsample of our clusters in the redshift range z = - 0.4. 

The observations made use of the EFOSC2 instrument at the 
3.6 m New Technology Telescope (NTT) in La Silla, Chile. The 
observations were carried out in September 2010, with typical 
exposure times of 840 seconds (two spectra per cluster, 420 sec- 
onds each). Our long-slit observations have been obtained us- 
ing Grism #4 (wavelength range 4085 - 7520A). The sHts (1.5" 
width) were placed on the BCG and an additional suitable cluster 
member candidate by rotating the slit. The BCG galaxies in these 
systems could be easily identified as the brightest red-sequence 
galaxies always coincident with the X-ray centroid, allowing us 
to safely anchor the cluster redshifts. 



A standard reduction process was applied to the data us- 
ing IRAF tasksQXhe observations were bias subtracted, cleaned 
from cosmic rays, and flat fielded. For each galaxy we have ob- 
tained two spectra which were sky subtracted and combined to 
increase the signal-to-noise ratio. The wavelength calibration 
was carried out by comparison with exposures of He and Ar 
lamps. 

The final spectra were then correlated with a database of 
galaxy templates. In almost all cases the H and K lines and the 
4000A break were visible and used for visual check of the tem- 
plate correction results. Spectroscopic redshifts have been se- 
cured for a total of 12 BCG galaxies. Due to relatively short ex- 
posures used only four systems a second member galaxy in the 
slit had good signal-to-noise ratio in order to safely measure its 
redshift. In all four cases the galaxies were found to have concor- 
dant redshifts with the BCG value. The spectroscopic redshifts 
of the galaxies are summarized in Table[3]along with our photo-z 
estimates. We compare the two redshift sets in Sect. 14.21 

4.2. Comparison of spectroscopic and pliotometric redsliifts 

For a subsample of 12 clusters (z < 0.4) we have obtained spec- 
troscopic redshifts of their BCG galaxy and in four systems also 
for one additional member galaxy (Sect. 14.1b . We compare the 
spectroscopic redshifts with our photo-z values in Fig. [8] (left). 
Our photometric values (red points) agree well within the error 
bars with the spectroscopic redshifts of the BCG (green points, 
brighter green points mark the clusters with two concordant red- 
shifts). Blue points mark the photo-zs for five of the systems 
obtained by the SCS survey (M09, MIO). These values exhibit 
a systematic bias toward lower redshifts, with a mean relative 
difference of 19%. A similar trend is also visible in Fig.[T3](top 
left) where we compare our photo-z values with the SCS mea- 
surements for all clusters common to both samples. 

The right panel of Fig. [8] displays the comparison of the ab- 
solute difference of our photometric and spectroscopic estimates 
in units of photo-z error, D = Izphoto -Zspecl/o"photo- A comparison 
with a Gaussian expectation shows an agreement at the 96% con- 
fidence level, confirming both the good precision of our photo-z 
estimates and the realistic description of their errors. 

The present spectroscopic sample covers only part of the red- 
shift range and does not allow us to check the photometric red- 
shift calibration at higher redshift. However, the good agreement 
at low z supports the photo-z method used. We also note, that 
same photo-z method has been applied to a large number of SPT 
selected clusters extending to beyond z = 1 , and the agreement 
between s pectro scopic and photometric redshifts has been excel- 
lent dHigh et allium Iwmiainson et aLliToTTl) . 

5. Results 

5. 1 . Galaxy cluster sample 

Table |2] provides the physical properties determined for the 46 
clusters in the present sample. The measured X-ray luminosity 
of the systems (Sect. I3.3l l and the photometric and spectroscopic 
redshifts (Sect. I4land l4. Il l are used as inputs for the cluster scal- 
ing relations to estimate further physical parameters. Ancillary 
X-ray information on the individual clusters can be found in 
TablelAl] 

The redshift, temperature and mass distributions are shown 
in Fig. |9] We display the X-ray luminosity of our systems as a 

iraf.noao.edu 
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Fig. 6. Left: Pseudo-color image of source ID 018 from the Blanco Cosmology Survey in the gri bands. X-ray contours are overlaid 
in white. The green circle shows the estimated 0.5 x rgoo. Right: Color-magnitude diagram {r - z vs. z) for the cluster ID 018. The 
red points show the member galaxies used for the photo-z estimate and the red dotted line indicates the single stellar population 
model at cluster redshift z - 0.39 (see text for details). 




Fig. 7. Pseudo-color images in the gri bands of the two X-ray detected (white contours) systems with secure photo-z values above 
z > 0.9. Green circles mark the 0.5 x rsoo radius. Both clusters have a BCG coincident with the center of the X-ray detection. 



function of redshift in Fig.[TO] The median redshift of the cluster 
sample is z = 0.47. Six of the systems have photometric redshifts 
z > 0.8. Three of these have redshifts consistent with z - 1, al- 
though the photo-z uncertainty in this regime is large. The me- 
dian temperature of the clusters is ~ 2 keV and the median M500 
mass 9 X 10'^ (based on luminosity scaling relations). We are 
thus able to probe the cluster/group transition regime practically 
at all redshifts out to z w 1 . 



5.2. Survey sky-coverage 

The simplest statistical characteristics of a cluster survey are its 
area coverage as a function of limiting flux {sky-coverage func- 
tion) and the cumulative surface density of the detected objects 
above the given flux limit as a function of flux - the so-called 
log - log S relation^ 



' We use the standard notation of this relation, but keep writing fx as 
the source flux rather than S . 
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Fig. 8. Consistency test of our photometric redshift estimates with spectroscopic measurements in the redshift range z - - 0.4. 
Left: Comparison of our photometric redshift estimates (red, Icr error bars) with spectroscopic values (green). Green stars mark 
clusters for which we have two concordant galaxy redshifts, while green circles indicate clusters for which only the BCG has a 
spectroscopic redshift. The photometric redshifts obtained by the SCS survey (M09, MIO) are shown in blue. The x-axis displays 
the cluster ID number. The objects are sorted in increasing redshift order The bottom panel shows the residuals of the photo-z 
values with respect to the spectroscopic measurement. Right: Cumulative histogram of the difference between the photometric and 
spectroscopic redshift normalized by the 1 cr uncertainty of the photo-z values, i.e. D - Izphoto - Zspecl/o"photo- The dashed line shows 
the expectation for the Gaussian distribution. Both curves are in good agreement, with a Kolmogorov-Smirnov test confirming that 
the distribution of the D values is Gaussian at the 96% confidence level. 



In order to properly determine the survey's sky coverage, 
good knowledge of the survey's selection function is necessary. 
For the simple case when the selection function is the function of 
only flux, the sky coverage is then the selection function of the 
survey scaled by its geometric area. Especially for the case of ex- 
tended sources the situation is more complex, since the selection 
function depends also on other parameters (e.g the source ex- 
tent and ofF-axis angle). These efifects can only be accounted for 
by Monte Carlo simulations. At this moment, without the sim- 
ulations at hand, we can still provide a preliminary, empirically 
calibrated sky coverage calculation and cluster log - log S re- 
lation. We will demonstrate that these simple approaches show 
good agreement with the design aims for the survey depth and 
previous measurements of the cluster log - log S function. 

While our source detection pipeline utilizes multiple detec- 
tion bands and likelihood thresholds (Sect. |3]l we will for sim- 
plicity (and ability to compare our results with published work) 
characterize detections made in the standard 0.5 - 2 keV band 
with a 3cr detection threshold and a 5cr extent significance. 

In order to obtain the survey sensitivity function for extended 
sources, we first calculate the point source sensitivity for each 
field. This is a simpler task since it does not require treatment 
of the source extent. We calculate the point source sensitivity 
function by analytically inverting the detection likelihood calcu- 
lation (described in Sect. 13.1b and obtaining the minimal count- 
rate necessary for a point source to be detected at the required 
detection threshold given the local background and exposure in 
the detection cell. 

The procedure is repeated for each survey field and the re- 
sults are combined for the whole survey area. In the areas where 
two or more fields overlap, we compare the sensitivity maps 
pixel-by-pixel taking the highest reached sensitivity (i.e. low- 
est local count-rate limit) at the given position. This procedure 
is chosen because the present catalog was derived from the de- 
tection pipeline that ran on each field individually. An alterna- 



tive approach is to combine the fields before detection - reaching 
slightly deeper flux-limits in the overlapping areas0This comes 
at the cost of losing the information on the local PSF shape used 
by the maximum-likelihood fitting algorithm, since the same sky 
location in two diff'erent observations is imaged at different off- 
axis and position angles and thus with different PSF. Both ap- 
proaches give comparable results and we opt here to characterize 
the main scheme (i.e. detection on individual fields). 

The median point source sensitivity calculated in this way 
for the whole survey area is 3.7 x 10"''' erg s"' cm"^ for an 
energy-conversion facto£| of 1.5 x 10"'^ erg s"' cm"^. The cor- 
responding sky coverage as a function of flux is displayed in 

Fig.m 

In the next step, we attempt to obtain a first order approxi- 
mation to the sky coverage function for the extended sources by 
a simple scaling to the point source function. In Fig. [10] (left) we 
show the dependence of the detection likelihood (i.e. the detjnl 
parameter) on the total detected source counts for point sources 
and the confirmed clusters from our sample (full circles). 

Photons from extended sources are distributed over a larger 
area and thus require more counts to reach a given detection like- 
lihood compared to point sources. For those, a simple linear re- 
lation in the log-log plane is a good description of the counts- 
detjnl relation (dashed red line in Fig. [TOb . Since the number 
of our clusters is small a similar linear relation for them is only 
very weakly constrained. We therefore fix the slope to the value 
from the point source fit leaving only the intercept as a free pa- 
rameter and weighting the points by their counts error (solid red 
line). The offset of the extended-source best-fit line translates to 



* This was done for the ancillary catalog using the wavelet detection 
algorithm. 

' Assuming a power law spectrum with Y = 1.7 and nu = 1.25 x 10^" 
cm"^ (median value of the galactic column density in the survey field) 
and using an on-axis PN response file. 
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Fig. 9. Top left: The redshift distribution of the 46 cluster sample based on the photometric redshifts ob tained with the red sequence 
fitting method. Top right: The X-ray temperature distribution estimated from the L-T scaling relation of'P ratt et al.l ( l2009l). Bottom: 
Distrib ution of the cluster masses in the rsoo aperture calculated from the luminosities using the L-M scaling relation from lPratt et al.l 
(|200§(see Sect. [331). 



a factor of 2.4 between the total required counts of point and ex- 
tended sources at any given detjnl. Fixing the slope has also 
the advantage that the offset is independent of the selected de- 
tection threshold. The best fit line roughly follows the locus of 
clusters with extent (beta model core radius) close to the median 
value of ~ 20". The solid red line in Fig. [10] (left) thus roughly 
gives the detected counts for a cluster with a typical extent de- 
tected with a given likelihood. We then use this off'set factor to 
scale up the point source sky coverage function (see Fig. fTTT i. 
The median flux limit for this sky-coverage is 9.3 x 10"'^ erg 
s"' cm"^ (using the median ECF of our sample). In Fig. [TOlwe 
display the luminosity-redshift plane for our survey. The lumi- 
nosity threshold for a flux limited sample (fmin = 1 x 10 erg 
s"' cm"^) is also shown, demonstrating a rough agreement with 
our calculation. Note that in the present sample we also include 
fainter sources than this threshold (the lowest cluster flux is 
~ 6 X 10"'^ erg s"' cm"^). 

This approach underestimates the effect of clusters with 
larger extent - and thus overestimates the sky coverage at given 
flux0 However, since the detection probability itself is a strong 



* The fit with a free slope gives an offset factor of ~ 4, the fit being 
skewed towards the locus of very extended clusters. 



function of source extent, the only way to properly account for 
its effect is through realistic simulations. 

We present examples of preliminary sky coverage functions 
for extended source detection on several (non-XMM-BCS) fields 
based on such Monte Carlo simulations (Miihlegger 2010) in 
Fig. IB. II discussed in Sect. |B] These first results validate our 
attempt to model to first approximation the extended source sky- 
coverage by scaling the point source curve and also confirm that 
the scaling factor between them is roughly ~ 2.4 (this scaling 
factor is expected to hold only for observations with roughly 
same depth as ours, ~ 10 ks). 

5.3. C/usfer log A^- log 5 

We now use this empirical sky-coverage in order to calculate the 
survey's log - log 5, defined in a standard way as: 



Nc 



N{> fx) = 



1 



deg 



(5) 



where Nc is the total number of clusters and f2(/^) is the ex- 
tended source sky-coverage corresponding to the flux of the i-th 
cluster. We characterized the survey sky coverage only for a hy- 
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Fig. 10. Left: Detection likelihood (detjnl) as a function of total detected source counts (PN detector only) for point sources (open 
circles) and the detected clusters (full circles). Clusters are color coded by their extent (beta model core radius). The dashed red line 
shows the best fit linear relation in the detjnl - counts plane for point sources. The solid line shows the same relation for extended 
sources (with slope fixed to the point-source fit, points weighted by their counts error). Typically, an extended source has to have 
2.4 times more counts than a point source to be detected at the same detjnl value. See Sect. l5.2l for details. Right: Luminosity in 
the 0.5 - 2 keV band (object rest-frame) for the present cluster sample as a function of redshift. The line shows the luminosity of a 
cluster with a measured flux of 1 x lO''* erg s"' cm"^ (unabsorbed, observer rest-frame). 




fx{0.5-2.0 keV) [erg s"'cm-^] 



Fig. 11. The survey sky coverage. Dashed line shows the sky 
coverage as a function of limiting point source sensitivity in the 
0.5 - 2.0 keV band. The empirically estimated extended source 
sensitivity is shown with a solid line. The median point source 
sensitivity of the survey is 3.7 x 10"'^ erg s"' cm"^, the median 
sensitivity for the extended sources 9.3 x 10"'^ erg s"' cm"^ 
(vertical line). 



pothetical single band (0.5 -2 keV) detection scheme. Since such 
a detection scheme is not part of our pipeline, we opt to draw a 
subsample from our cluster catalog derived from the three band 
scheme (which includes the 0.5 - 2 ke V). For this calculation we 
consider only clusters that would have also been detected in this 
hypothetical single band run by setting the same detection- and 



extent likelihood thresholds used for the sky-coverage calcula- 
tion in the previous section. 

This requires us to recover the actual 0.5 - 2 keV band detec- 
tion likelihoods from the total detjnl parameter, which includes 
contributions from all three detection bands. As we described in 
Sect. |3] detjnl can be interpreted as detjnl = -InPrand, with 
Piand being the probability of a false detection arising from pure 
Poissonian fluctuations. The actual definition of this parameter 
is slightly more complex: 

detjnl = - ln(l - r(0.5v, L)), (6) 

where F is the incomplete gamma function and its arguments 
are the number of degrees of freedom of the emldetect fit (for 
extended sources v = 3+ the number of detection bands times 
number of instruments) and L is t he sum of a ll the individual 
likelihoods (using the C statistics ofl CashlFl979l) . This definition 
effectively converts the joint likelihoods to two degrees of free- 
dom allowing to compare detections from different combination 
of bands and instruments. However, for the conversion to a sin- 
gle band detection likelihood, we need the original individual 
likelihoods which we obtain by numerically inverting Eq. |6]for 
each source. 

We then calculate the joint detection likelihood from all three 
instruments in the single, 0.5 - 2 keV, band (all three detec- 
tion probabilities being independent) and subsequently calculate 
the new detjnl parameter normalized back to two degrees of 
freedom using Eq. |6] The number of clusters that have this new 
single band detjnl parameter above the required threshold (i.e. 
equivalent to ~ 3cr) is 40. 

Finally, we calculate the logA^ - log 5 according to Eq. |5] 
and the variance of the number counts as cP- - Yl^fi l/Q(/p^. 
The recovered curve (see Fig.fT2ll is in good a greement with the 
logAf - log S of other surveys: e.g. COSMOS (iFinoguenov et al.l 
l20()7h . the RDCS survey (the ROSAT Deep Cluster Survey, 
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Fig. 12. The log - log S of the present sample in the 
0.5 - 2.0 keV band. Fluxes are calculated in the rsoo aper- 
ture. Res ults from several surveys are also shown: COSMOS 
dFinogue nov et al. 2007), RDCS (the ROSAT Deep Cluster 
Surve y, Rosati et al. 1998), 400 deg^ survey ("Bure nin et al 
2007t " rVikhlinin et al. 2009) and the XMM-LSS (iPacaud et al 



2003). See Sect.lSJlfor details. 



Rosati et al.' '1998"), 400 deg^ survey (iBurenin et al.' '2007^, 
Vikhlinin et al. 2009) and the XMM-LSS surve>0 (Pacaud et al. 
2007h . Since the area and depth of the XMM-LSS survey match 
well the parameters of our survey we discuss their comparison 
in more detail in Sec. IC.ll 

We note, that we used the fluxes in the rsoo aperture for our 
calculation whereas the XMM-LSS uses a fixed physical aper- 
ture of 0.5 Mpc (typically very close to rsoo), RDCS ~ 80 - 90% 
of the total flux (i.e. integrated out to infinity) and the 400 deg^ 
survey the full total flux. We have chosen rsoo because: a) it re- 
quires less extrapolation based on a beta model whose param- 
eters are typically highly uncertain and is itself not necessarily 
a good description of the surface brightness profiles and b) it is 
the most natural choice when comparing to theoretical predic- 
tions (the cluster is approximately virialized in this radius and 
the scaling relations we employ are calibrated for this overden- 
sity). However, assuming a typical cluster with rjoo - 0.5 Mpc 
well described by a beta model with (J3, rcore) - (2/3, 180 kpc) 
the flux extrapolated to infinity would be higher by ~ 1/3 mov- 
ing our curve along the x-axis to higher fluxes only very slightly 
- even closer to the RDCS and 400 deg- survey's relations. 

Uncertainties of the flux estimation (including the uncer- 
tainty on the photo-zs) affect the logA^ - log 5 only in a mi- 
nor way. The main source of uncertainty (not included in the 
error bars) is our current lack of knowledge of the survey se- 
lection function (and thus only tentative description of the sky- 
coverage). The good agreement with previous work gives, how- 
ever, support to our preliminary approach. 

5.4. Cross-correlation with known sources 

The XMM-BCS field has an extensive multi-wavelength cover- 
age and has already been studied by the Southern Cosmology 



Survey (M09, MIO. iMcInnes et al.ll2009l) who identified in op- 
tical data a number of clusters in this area. Due to a significant 
overlap with our cluster catalog we will address a more detailed 
comparison in Sect. 15.5! 

In search of other known sources coincident with our 
clusters, we make use of both the NASA/IPAC Extragalactic 
DatabasfBand the SIMBAD Astronomical Database^ 

First we looked for associated known clusters. For this 
query a search radius of 60" was selected, finding a sin- 
gle match - the cluster 400d J2325-5443 (alternative name: 
[BVH2007] 240) identified in the 160 Square De gree ROSAT 
Survey (Vik hlinin etaP 119981: iMulhs et al] l2003h at spectro- 
scopic redshift z - 0.102. This cluster is coincident with our 
cluster ID 476 with photometric redshift of 0.1 being in full 
agreement with the spectroscopic value. The source is also part 
of the 400 Square Degree ROSAT Survey. See Appendix lA!2l for 
more details on this source. 

We also list galaxy matches, if they are within a 16 arc- 
sec search radius from the X-ray center in Table IA.2I (in the 
appendix) with matches coming from the 2 Micron All Sky 
Survey Extended objects catalog and the APM galaxy survey 
(ISkrutskie et al.ll2006LlMaddox et alJI 19901 respectively). Out of 
13 matches, only two galaxies have known spectroscopic red- 
shifts, both obtained in the frame of the 6dF Galaxy Survey 
(IJones et al.lf2004l) . The first is 2MASX J23254015-5444308 at 
z = 0.101 coincident with the brightest galaxy in cluster ID 
476. The redshift value is concordant with the redshift from 
the 160/400 Square Degree ROSAT surveys. The second match 
is the brightest cluster galaxy of the system ID 150 at redshift 
z - 0.176, again in good agreement with our estimated photo-z 
of 0.2. 

As can be seen, the survey field has a wealth of multi- 
wavelength data, but very little spectroscopic measurements. 
This makes the ongoing spectroscopic follow-up program very 
important, as redshifts are essential for the full utilization of the 
available data sets. 

Radio sources coincident with the X-ray detected clusters 
can bias the SZE signal (filling the decrement). We checked for 
intervening radio sources by cross-correlating our cluster catalog 
with the NED database with a 1 arcmin search radius. We find 
1 1 radio sources detected at 843 MHz by the Sydney University 
Molonglo Sky Survey (SUMS, MauchetaL 2003). The source 
PMN J2330-5436, ~ 30 arcsec from cluster ID 150, was detected 
by the Parkes-MIT-NRAO (PMN) southern survey at 4.85 GHz 
(Gregorv et al. 1 994) . The list of all identified radio sources is 
given in Tabled 



5.5. Cross-correiation witti ttie Southern Cosmoiogy Survey 
dusters 

The Southern Cosmology Survey (SCS) carried out an optical 
cluster search using the Blanco Cosmology Survey imaging data. 
Menanteau et al. ( 2009, hereafter M09) provided a catalog of op- 
tically selected clusters with photo-z < 0.8 in a 8 deg^ field par- 
tially overlapping with the 6 deg^ region presented in this work. 
iMenanteau et al. (1201 Oi MIO) then followed up this work by cre- 
ating a cluster catalog using the full 2005 - 2007 BCS survey 
data (i.e. ~ 70 deg^, thus fully covering also the whole XMM- 
BCS field), detecting 105 clusters with M200 > 3 x lO'"* Mq and 
photo-z < 0.8. 



^ Note that the XM M-LSS curve is only digitized from the figure 
in IPacaud et al.l d2007l) since the original curve is no longer available 
(Pacaud, private com.). 



nedwww.ipac.caltech.edu 
" simbad.u-strasbg.fr/simbad/ 
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Fig. 13. Top left: Compariso n of photometric redshifts for the 19 common clusters from our sample and the SCS cluster survey 
from iMenanteau et alj (l2009l) and iMenanteau et al] (1201 0). Top right: Comparison of masses for the same cluster sample in the 
r2oo aperture determined from the measured X-ray luminosity through scaling relations (x-axis) and the mean of the optically 
determined masses M(L2qq) and M(N2oo) (Table|5]l. The red line marks equality in both top panels. Bottom left: Photo-z difference 
A = photo-z(SCS) - photo-z(XMM-BCS) as a function of our estimates of redshifts. Bottom right: Mioo difference A = M2oo(SCS) 
- M2oo(XMM-BCS) as a function of our estimates of redshifts. The optical masses are significantly higher than the X-ray estimates 
especially at the low and high redshift ends. See text for discussion. Green points in all plots mark clusters from MIO, black points 
those from M09. 



Combining both these catalogs, we find in total 30 SCS clus- 
ters whose optical coordinates lie inside our 6 deg^ regionQOut 
of these 30 systems, 26 come from the M09 catalog (which con- 
tains clusters also below the mass limit applied in MIO), two are 
included in both M09 and MIO and two additional clusters are 
from MIO. For the two clusters in both M09 and MIO we will 
use the updated parameters from MIO. 

The SCS catalog provides the BCG coordinates for each 
system while our catalog lists the X-ray centroids. For cross- 
correlation of the two catalogs we take a conservative 60" 
matching radius, which yields 19 clusters. We summarize the 
properties of matched clusters in Table |5] and provide a more 
detailed comparison of their parameters in the next two sections. 



In the present work we thus do not consider SCS clusters that lie 
only partially in our 6 deg^ region or in the 8 deg^ mosaic extension. 



5.5.1. Comparison of photometric redsliifts 

First, we turn to the comparison of the photometric reds hifts for 
the 19 matches. M09 and MIO utilize the BPZ code dBemtezI 
2000) to estimate photo-zs while our method is based on the red- 
sequence method as described in Se ct.lH For the SCS clus ters we 
use the photo-z errors published in 'M cInnes et alj (l2009l) where 
possible (the M09 and MIO catalogs do not provide error bars). 
For the remaining cases we assume a 15% error, which is the 
mean precision of the photo-zs where errors are available. 

As can be seen in Fig.[T3](top left) there is no case of catas- 
trophic disagreement. We find a gap in the SCS photo-z distribu- 
tion in the 0.35 - 0.5 photo-z range that is not present in our 
redshift distribution. The most important feature is, however, 
the systematic offset between the photo-z estimates. The SCS 
photo-zs are on average ~ 20% lower than our values. This trend 
roughly holds in the whole redshift range, as can be seen from 
the photo-z residuals plotted against redshift in Fig. [T3] bottom 
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Table 4. Radio sources within 60" from the X-ray centers of the detected clusters. The quoted flux density S is at 843 Mhz (36 
cm) for the SUMSS sources and at 4.85 GHz (6.2 cm) for the PMN detected object. The radio counterparts were obtained from the 
NASA/IPAC Extragalactic Database. 



ID 


Object Name 


R.A. (deg) 


Uec (deg) 


S (mJy) 


separation (") 


018 


SUMSS J232952-560723 


352.4677 


-56.1231 


1 /I n 

14. y 


± U.8 


56 




QTTA/TQQ T^^^^A^ ^'>^S17 


QCQ AA^ A 


-J J.UJOZ 


16.7 


±0.8 





044 


SUMSS J23 1654-545406 


349.2274 


-54.9019 


8.0: 


t 1.0 


14 


109 


SUMSS J232737-541622 


351.9047 


-54.2730 


26.1 


± 1.2 


9 


110 


SUMSS J233003-541424 


352.5146 


-54.2402 


13.3 


± 1.1 


6 


189 


SUMSS J233044-560123 


352.6860 


-56.0233 


15.1 


±0.8 


36 


210 


SUMSS J233406-554708 


353.5253 


-55.7857 


7.9: 


tO.7 


3 


288 


SUMSS J233459-545535 


353.7495 


-54.9265 


41.5 


± 1.6 


37 


426 


SUMSS J232138-541849 


350.4092 


-54.3137 


12.7 


± 1.9 


14 


534 


SUMSS J232446-552432 


351.1951 


-55.4089 


41.1 


± 1.5 


17 


546 


SUMSS J233 113-543025 


352.8076 


-54.5071 


25.6 


± 1.3 


26 


150 


PMN J2330-5436 


352.5075 


-54.6097 


52.0 


±8.0 


34 



left). We found a similar bias when comparing the SCS photo-zs 
to the spectroscopic subsample in Sect. 14.21 For the five systems 
with spectroscopic redshifts the photo-zs were on average un- 
derestimated by ~ 19%. 

In order to investigate potential sources of the discrepancy 
we check whether the photo-z residuals depend on any of the 
available parameters, most importantly the richness parameter 
A^200, integrated optical luminosity and the BCG-X-ray cen- 
troid offset. However, we do not find any statistically significant 
dependence. 

5.5.2. X-ray - optical mass comparison 

M09 and MIO provide rough mass estimates for their clusters 
based on t he optical proxies N ^m and Ljg^ using the scaling rela- 
tions from'Reves et al.l (1 2008 ). The richness estimator Af2oo is de- 
fined as the number of E/SO ridgeline cluster members brighter 
than 0.4L*. The integrated cluster luminosity L°pg is the summed 
r band luminosity of the member galaxies included in the A^200 
calculation. Both parameters are calculated within an aperture 
where the galaxy density equals 200/Qm t imes the mean density 
of galaxies in the Universe. Fortunately, Johnston et al.l (l2007h 
found that this aperture is an unbiased estimate of the radius 
where the matter density is 200 times the critical density of the 
Universe, i.e. the optical masses and our X-ray estimates come 
from roughly the same apertures and can thus be directly com- 
pared. 

In Fig. [13] (top right) we compare our X-ray masses with 
the optical masses M(A^20o) calculated from the A^200 parameter. 
The optical masses are estimated to be accurate within a factor 
of two (M09), where this f actor should includ e also the uncer- 
tainty in extrapolating the Reves et al] (|2008|) scaling relations 
to higher redshifts (the scaling relations were calibrated for red- 
shifts z < 0.3). We used the factor two uncertainty to calculate 
the M(A^20o) error bars in Fig. [13] We find that the optical masses 
are significantly higher than the X-ray mass estimates M200 by a 
factor of ~ 2.6 (median value). 

iReichert et al. (2011) investigate X-ray luminosity based 
scaling relations on a large compilation of cluster samples from 
the literature. They find only very few systems deviating from 
the mean L - M relation by more than a factor two (i.e. with 
actual mass two or more times higher than the luminosity pre- 
diction). We thus do not expect our X-ray masses to be under- 
estimated by similar factors even in individual cases. The ob- 
served bias in mass goes in the opposite direction as that found 



in the photo-zs (i.e. photo-zs were underestimated while masses 
overestimated). The photo-zs, however influence the mass esti- 
mates and therefore it is not straightforward to disentangle all 
the factors contributing to this discrepancy. The influence of the 
redshift uncertainty is likely more important for nearby systems, 
where it translates to larger differences in the angular size of the 
aperture. The discrepancy for the M{L^^^ masses is similar. We 
note here, however, that the M{L°^ ^ masses iri M09 were ob- 
tained from the scaling relations of lReves et al.l (l2008h prior to 
their erratum-correction[3 

The bottom right panel of Fig. [T3]displays the mass residuals 
versus our photometric redshifts. We find an anticorrelation be- 
tween redshift and the mass residuals (using a Spearman's rank 
correlation c oefficient, neglecting the error bars). Given that the 
iReves etal] (j2008) relations are calibrated only out to z < 0.3 it 
is however impossible to either confirm the presence of such an 
anticorrelation or its potential causes. We also note that the four 
clusters from MIO agree with our measurements better than all 
but one cluster from M09 . 

iMcInnes et alJ (l2009h provides weak lensing mass measure- 
ments for the clusters from M09 among which 13 clusters are 
also in our sample (for two of these systems only upper limits 
could be set). We compare the weak lensing masses with our 
X-ray estimates in Fig. [14] The agreement is significantly better 
than for M(A^20o) masses, although the scatter and uncertainty 
in the weak lensing mass estimates are large (mostly due to the 
limited depth of the optical data, but also influenced by the un- 
certainty in the photo-zs). From their full sample. iMcInnes et al.l 
^09) also noted that the M(L°p|,) seems to overestimate the to- 
tal mass compared to the weak lensing estimates. 

An in-depth comparison of optical and X-ray masses will be 
addressed in an upcoming work, where we will provide also our 

own measurements of N200 and l^^m ('^o'^g ^1-' prep.). This 
will allow us to properly investigate the presence of potential 
biases in the different mass estimators methods and cahbrate our 
own relations. 

5.5.3. Parameter upper limits for X-ray non-detections 

For the 1 1 SCS clusters that lie in the core area of our survey but 
have no X-ray counterparts we provide X-ray flux upper limits 
in Table|6]and mass limits using the SCS photo-z values. 



Scaling relations with the updated coefficients are available at: 
arxiv.org/abs/O802 .2365. 
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Table 5. The 19 SCS clusters recovered in the XMM-BCS survey. SCS References: M09 - iMenanteau et al.l (l2009h . MIO - 
[Menanteau et al.' ("20101). The optical mas ses M(L'I^^q) and M(N2oo) are taken from M09 and MIO, the weak lensing mass mea- 
surements are provided bv lMcInnes et al.l (2009). The X-ray mass estimates obtained in the present work are taken from Table |2] 
The ID of the X-ray counterpart and its distance from the BCG are listed in the last two columns. The masses are in units of 10'^ Mq. 
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Table 6. The 1 1 SCS clusters from IMenanteau et alJ (l2009l) that lie in the XMM-BCS core survey but have no X-ray cluster 
detection, for which we provide flux and mass upper limits. 
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The flux limits were calculated using the same procedure as 
we used for the survey sky coverage calculation (Sect. 15.21 i.e. 
calculating the minimal flux needed for the source to be detected 
at the given position and our detection threshold). 

The flux was then converted to luminosity using the pho- 
tometric redshift from either M09 or MIO. We calculated the 
mass upper limits from the L - M scaling relation as detailed in 
Sect. l3.3l The obtained upper limits on the mass are considerably 
lower than the M09/M10 estimates. 

We also check for possible miss-classification (or confusion 
if a central AGN is present) by cross-correlating the positions 
of these 1 1 clusters with our X-ray point source catalog with a 
threshold of 16". In this aperture we find no matches. The non- 
detection of these sources is due to either thek low X-ray fluxes 
or a spurious detection. 



6. Discussion 

In this section we discuss the additional effects that influence the 
precision of the physical parameters provided in our catalog. We 



also give an outlook on the upcoming work in the context of the 
XMM-BCS survey. 

6.1. Error budget of the X-ray analysis 

For the present catalog, we restricted ourselves to provide only 
formal statistical errors for the estimated parameters (Table |2]) 
that include the Poisson errors of the flux measurement, a 5% 
systematic error from the background modeling and the intrinsic 
scatter of the scaling relations. Although we used the bolomet- 
ric luminosity to calculate further physical parameters, here we 
assumed the intrinsic scatter found in the 0.5 - 2 keV luminos- 
ity relations. This scatter is slightly larger than the bolometric 
one and it gives a more realistic error estimate since the band 
luminosity is, in fact, our only direct observable, while the tem- 
perature required for the bolometric correction is not. We de- 
termine physical parameters with the following precision (mean 
across the whole redshift and flux ranges): flux and luminosity 
to ~ 16%, T500 and Mgoo to ~ 30%, and Ygoo to ~ 60%. 

In this section we discuss several additional sources of sys- 
tematic errors and their impact on the estimated fluxes and other 
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Fig. 14. Comparison of X-ray masses (M2()o, x-axis) with 
the weak le n sing measurements (M(WL), y-axis) from 
iMcInnes et al.l (12009^ for 13 clusters in our sample. Although 
the scatter and uncertainties are large, the agreement is consid- 
erably better than with the optical masses M(A^20o) displayed in 

Fig.ini 



parameters. All below reported relative errors are obtained by 
averaging over the whole cluster sample. Several of the consid- 
ered effects are redshift dependent, but we typically allow broad 
parameter ranges and thus our uncertainty estimates are rather 
conservative. 

1) Good precision photometric redshifts are crucial for the 
determination of each physical parameter Photo-z estimates 
in the present work have a mean error of ~ 10% and show 
good agreement with the available spectroscopic measurements 
(Sect. 14.21 ). In order to estimate the impact of the photo-z un- 
certainty on the measured physical parameters we offset the red- 
shifts (Table |2]l by their Icr errors to both sides and rerun the it- 
erative physical parameter estimation procedures (see Sect. 13.31 ). 

We find that, for the flux fx(< /"500), all values are consis- 
tent within their 1 cr uncertainty and for most clusters the rela- 
tive difference is below the ~ 2% level (Fig. [15] left). Change 
in the photo-z affects the flux in a complex way - it is entering 
directly the energy-conversion factor (ECF) calculation (lower 
redshift leads to a lower ECF), and also during the iterative pro- 
cess through the scaling relations, which then feed back into the 
aperture size itself as well as the temperature which again affects 
the ECF value. This complex dependence explains the scatter of 
the flux residuals in Fig. [TSj leading to different convergence 
points for different input photo-zs. Interestingly, a lower photo-z 
value leads to a higher flux in the rsoo aperture. The reason is that 
the direct effect of decreasing the photo-z would be to lower tem- 
perature and mass and thus also reduce the rsoo value. However, 
the redshift dependence of the angular distance is stronger and 
thus the angular size of the rsoo aperture is actually larger for 
lower redshifts, which leads to the increase in the fx(< rsoo) val- 



ues (we confirm this explanation by checking the flux in fixed 
sky apertures). 

For luminosities the photo-z errors translate into a ~ 20% 
uncertainty (Fig. [15] right). Here the dependence is dominated 
by the cosmological redshift dimming and thus higher redshifts 
yield also higher luminosities. If we now use the perturbed 
redshift and luminosity values to recalculate temperatures and 
masses, we find that the T^qq values vary on the ~ 5% level, 
while for M500 the uncertainty is on the ~ 7% level. 

2) In the present work, we have utilized the bolometric lu- 
minosity scaling relations of Pratt et al. (2009) based on the 
REXCESS cluster sample (Bohringe r et alj 120071) . These scal- 
ing relations are best suited for our purposes for several reasons. 
They were derived from XMM-Newton observations (removing 
possible calibration issues between relations derived from differ- 
ent instruments) of a representative cluster sample (REXCESS, 



Bo hringer 



gents) 
et alJ [ 



2007). The sample covers a great range of the 



cluster luminosity function without a bias towards a morpholog- 
ical structure type (like e.g. presence of a cooling core or merg- 
ing activity). By using the bolometric relations we can also uti- 
lize the results of Reichert et al. (2011) to estimate the effect of 
non-self-similar evolution on the estimated parameters (see be- 
low). Additionally, the L- M relation is based on the L - Fx and 
Yx - M scaling relation, which is found to be more r obust than 
previous direct L - M calibrations (Arnau d et alj|200 7). 

The direct application of these scaling relations, however re- 
quires extrapolations both to higher redshifts (the REXCESS 
cluster sample includes only local clusters with z < 0.2) and to 
the low-mass regime of groups of galaxies. The physical param- 
eters provided in Table|2]were obtained by assuming the redshift 
evolution of the cluster scaling relations to be self-similar. This is 
a standard assumption supported by pre dictions of a p urely grav- 
itationally driven cluster growth (e.g. lKaise^ll986h . However, 
there is increasing evidence, that the evolution of the luminos- 
ity s caling relation s is s lower than the self-similar expectation 
(see iReichert et alJ (1201 li) and references therein). We test the 
influence of this assumption by usin g the simpUfied approach 
proposed bv'Fassb ender et al] (1201 Ibl) by removing a part of the 
self-similar evolution factor from the relations (the L - T "no 
evolution" scenario). This approach is consistent with the more 
detailed analysis of Reichert et al. In this picture, the predicted 
temperatures are on average higher by 9% and masses by 15% 
compared to the self-similar scaling relations. At the high red- 
shift end (z > 0.8) this effect is even more important (~ 20% and 
30% increase, respectively), while for z < 0.2 the effect is less 
than 5%. 

The cluster sample of IVikhlinin et alJ (l2009l) covers a sim- 
ilar redshift range as our sample (0 - 0.9) and extends in mass 
down to roughly the median mass of our sample (~ 1 x 10'** Mq). 
Using the L - M relation from this samples gives masses only 
slightly higher (by ~ 5%, standard deviation ~ 13%) than our re- 
sults. This difference grows with redshift due to the deviation of 
the evolution index from the self-similar value dVikhUnin et al.l 
(2009) obtained ^(z)'-^'). 

Pertaining the assumptions that the REXCESS s caling re- 
lations hold also beflow M500 ~ 1 x lO''* Mq, ISun et all 
(2011) found very good agreement between the pressure 
profiles of their sample of groups and those derived from 
REXCESS. i Leauthaud et all (1201 Q) used the gro ups detected 
in the COSMOS survey dFinoguenov et alj l2007h to caHbrate 
the L - M scaling relation with weak lensing mass measure- 
m ents. Their evo l ution factor is consistent with the one found 
bv lReichert et alj (1201 lb . i.e. giving only slightly higher masses 
than the self-similar scenario. These findings give us an indica- 
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Fig. 15. Left: The effect of the photometric redshift uncertainty on the determined flux in the rsoo aperture. Right: The effect of the 
photometric redshift uncertainty on the determined luminosity in the rsoo aperture. Green points mark the relative difference of flux 
(luminosity) for photo-z increased by 1 cr compared to the mean value. Red points are for the case when we decrease the photo-zs 
by the same amount. 



tion that the scaling relations employed in the present work are 
indeed adequate also in the low-mass clusters/groups regime. 

3) The Galactic hydrog en column densities reported by the 
LAB HI survey (iKalberla et al. 2005h are system atically lower 
by ~ 27% than the Dick ev & Lockmanl ( 119901) values in the 
whole survey area. The effect on the derived luminosities is, 
however, only marginal (~ 1.6%). 

4) In order to quantify the effect of possible deviations of 
the cluster metallicity from the mean value of 0.3 solar, we 
bracket the possible metallicities in the very conservative range 
of (0.1, 0.6) solar The corresponding range of Lx(< '"soo) de- 
viations from the fiducial value (for 0.3 solar metallicity) is 
(1.6%, -1.5%), i.e. lower metallicities lead to higher luminosi- 
ties and vice versa. 

5) We also test the quality of the flux extrapolation correction 
described in Sect. 13.31 The correction coefficients are calculated 
by integrating a beta model between (r^/o,, rsoo), if rpiat < r^oo. 
We use Eq.[T|to estimate the /3 and rc parameters. Alternatively, 
we can use the canonical value yS = 2/3 and rc obtained from 
the maximum likelihood fit in the source detection step. The two 
extrapolation methods give fluxes (luminosities) differing on av- 
erage for the whole sample by 2%. For individual objects the 
relative difference of fluxes is clearly correlated with the amount 
of extrapolation needed and is roughly of the size of the correc- 
tion itself. This means that the extrapolation is currently only 
very weakly constrained and thus highly uncertain. Fortunately, 
for most of our sample the required extrapolation factor is small. 

6) The combined MOSl and M0S2 counts are converted 
to flux and luminosities using the MOS2 response matrix (see 
Sect. 13.3b . We have chosen the M0S2 response matrix over the 
MOSl, because some sources lie on the missing MOSl CCD#6, 
where no meaningful response matrix can be calculated. If the 
MOS 1 response is used instead, the luminosities obtained purely 
from the combined MOS detectors are on average lower by 
2% (excluding clusters detected on the position of the missing 
MOSl chip). The final Lx(< '"500) calculated as the weighted av- 
erage of the individual PN and MOS luminosities is affected by 
less than 1.3%. 



7) The response matrices used in our analysis are calculated 
for a fixed radius of 150 arcsec. This range is roughly the av- 
erage extraction radius of our clusters (i.e. from which growth 
curves are extracted and local background estimated). We cal- 
culate response matrices for two additional radii - 60 arcsec and 
240 arcsec, to check how the spatial averaging of the spectral re- 
sponse impacts the derived ECFs and thus flux and luminosity. 
In this very conservative range of extraction radii we found the 
average effect to be of the order of 2.5%. 

8) The uncertainties in the absolute normalization of the 
effective area of the detectors decrease the flux measurement 
precision. INevalainen et al.l (i2010) found an agreement between 
0.5 - 2 keV fluxes measured by PN and both MOS cameras to 
be better than ~ 5 - 7% and for ACIS on Chandra and the PN 
found the fluxes to differ only by 2%. 

9) We have tested the impact of using the up-to-date version 
of the SAS package (SAS v. 11.0.0), with respect to the SAS 
V. 7.1.0 adopted for the analysis in this paper In an end-to-end 
reprocessing of a subsample of six clusters fully sampling our 
flux and off-axis angle distributions we find on average only very 
small differences (< 4% in flux, < 1.5% in temperatures and < 
2.5% in masses). These differences can get larger for low quality 
candidates (e.g. quality flag 3 in Table lA.ll . but in all test cases 
they were consistent within the error bars. 

10) We also tried to run the physical parameter estimation 
procedure (Sect. l3.3T l from several initial values of T500/ rsoo. The 
iteration procedure always converged to the same solution, con- 
firming its independence from the starting values. 

We will provide tests of photometric accuracy of the growth 
curve method in a subsequent publication based on simula- 
tions using realistic backgrounds (i.e. using our survey fields as 
background for the simulated clusters). Presently we provide a 
comparison of the X-ray photometry obtained by our pipeline 
in comparison with the XMM-LSS project in Appendix IC.ll 
In addition, our algorithm was also applied to the clusters 
SPT-CL J2332-5358 and SPT-CL J2342-5411 in ISuhada et aH 
(2_Q10) detected in the extension of the XMM-B CS survey. These 
sources have been independently analysed by lAndersson et al.l 
(1201 Oh using deeper pointed data (from a 19 ks long XMM- 
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Newton observation for the first cl uster and from a 1 34 ks long 
Chandra exposure for the second'). [Andersson et al ] (12010^ find 
for SPT-CL J2332-5358 (Zphoto = 0.32) an X-ray luminosity 
Lgoo = 3.0 + 0.1 X 10"^ erg s"' and for SPT-CL 2342-5411 
(Zphoto - 1 -08) Lgoo = 2.9 + 0.3 X 10"^"* erg s"' in good agreement 
with our values: L500 = 2.7+0.2x10'^'* erg s"' for SPT-CL 2332- 
5358 and Lsoo = 2.9 + 0.3 x lO'*'* erg s"' for SPT-CL 2342-541 1 
(luminosities in the 0.5 - 2 keV band). Estimated temperatures 
and masses are also consistent within error bars although the un- 
certainties on these parameters are significant given the exposure 
time in the Suhada et al. (2010) analysis was < 3 ks. 

A proper understanding of a realistic error budget of a cluster 
sample is crucial for its modelling in the cosmological context. 
From our analysis we find that most effects are typically on the 
~ 2% level (under conservative assumptions) and the major con- 
tributing factors are the uncertainty of the photo-z measurements 
and the required extrapolations of the scaling relations (both in 
the range of 5 - 30% depending on the parameter and the red- 
shift of the system). For a few clusters an additional significant 
source of uncertainty is connected with the flux extrapolation. 
A full self-consistent treatment of the error propagation (includ- 
ing their full covariance matrices) and its impact on the cosmo- 
logical modeling of the sample will be addressed in subsequent 
work. 

6.2. Project outlook 

The present sample establishes the observational base of the X- 
ray part of the XMM-BCS survey. In upcoming work we will 
use the available multi-wavelength data to follow several lines 
of investigations, some of which have already been initiated: 

- The X-ray cluster catalog will be extended to cover the whole 
14 deg^ area. The preliminary source catalog is already avail- 
able and we will follow up this work by estimating the pho- 
tometric redshifts and physical parameters for the clusters in 
the same way as presented in this work. The full cluster cat- 
alog is expected to comprise ~ 100 clusters and groups of 
galaxies. 

- We will calculate the selection function bas e d on Monte 
Carlo simulations developed by iMiihlegsie^ (l2010l) . This 
analysis will allow us to construct a well controlled subsam- 
ple from the full cluster catalog that will be suitable for cos- 
mological modelling. 

- A more detailed analysis of optical properties of the clusters 
presented in this sample wiU be provided in Song et al., in 
prep. We will provide here measurements of the A^200 and 
Ljpg parameters and investigate their mass scaling relations. 

- A detailed comparison of the X-ray, optical and mid-infrared 
cluster samples will allow us to gain good understanding of 
the selection function of each method. We will study the 
cluster/group population in this field and establish its multi- 
wavelength properties. The Spitzer imaging data will also 
be used to improve the photometric redshift estimates, espe- 
cially for distant systems with redshift z > 0.8. 

- We have initiated further X-ray-SZE studies based on a co- 
operation with the SPT co llaboration. The c urrent SPT clus- 
ter sample jVanderlinde "et al. 201 0: Wi lliamson et al.ll2011l) 
includes only sources with minimal detection significance 
of 5cr (7cr), respectively. There are only two clusters in the 
14 deg2 above this threshold, SPT-CL J2332-5358 and SPT- 
CL J2342-5411, which are also independently detected in 
our survey (iSuhada et alJl2010l) . Using our X-ray selected 



cluster catalog we can also safely investigate lower signif- 
icance SPT detections. As a first example, cluster ID 044, 
(XBCS 231653.1-545413) was found to have a direct SPT 
detection at the 4.2cr level (B. Benson, private communica- 
tion). Another approach is a stacking analysis of the SZE 
data for the X-ray selected clusters. Here, a preliminary anal- 
ysis of the top eleven clusters ranked by their X-ray pre- 
dicted SPT detection significance yields a > 6cr detection. 
We will explore both approaches in more depth in upcoming 
work, but already now it is clear, that with a joint SZE and 
X-ray analysis we are able to explore a completely new mass 
regime in the SPT survey. 

- The multi-wavelength coverage of the field provides oppor- 
tunities also for non-cluster science. As an example, we have 
detected a total of 3065 X-ray point sources in the survey 
(1639 in the core region and 1426 in the extension). Most of 
these point sources are AGN and using the available multi- 
wavelength data we will be able to carry out a study with a 
focus on the obscured AGN population. 

7. Summary and conclusions 

- We have provided the analysis of the 6 deg^ XMM-A^ewfo« 
field in the framework of the XMM-BCS survey. We have 
carried out X-ray source detection and constructed a catalog 
of 46 clusters and groups of galaxies. 

- Based on four band optical imaging provided by the Blanco 
Cosmology Survey we have confirmed that these X-ray de- 
tections are coincident with overdensities of red galaxies. 
Using the red sequence method we have measured the pho- 
tometric redshifts of these systems. 

- We have initiated a spectroscopic follow-up program by 
carrying out long-slit spectroscopy observations using the 
EFOSC2 instrument at the 3.6 m NTT telescope at La Silla, 
Chile. We have obtained spectroscopic redshifts for BCG 
galaxies in 12 clusters and in four cases also for one addi- 
tional member galaxy. This sample covers the redshift range 
< z < 0.4 (i.e. roughly up to the median redshift of the 
sample) and constitutes the first spectroscopic information 
for the field. We find good agreement between our photo- 
metric estimates and the spectroscopic values, but the spec- 
troscopic sample has to be extended in redshift, in order to 
be able to provide a rigorous calibration of the photo-zs. 

- Using the redshift information we measured the X-ray lumi- 
nosities for our cluster sample. From luminosity scaling rela- 
tions we estimate their most important physical parameters, 
e.g. mass, temperature and the Yx parameter We discuss the 
influence of several factors on the precision of the provided 
estimates. The uncertainty of the photometric redshift esti- 
mates and the extrapolation of the scaling relations to high 
redshift systems and into the group regime are identified as 
the most important factors that determine the overall errors in 
the physical parameters. We verify our X-ray parameter es- 
timation met hod by analyzing th e CI sample of the XMM- 
LSS survey (iPacaud et alJl2007h . We find good agreement 
between the parameters provided by both pipelines. 

- The present sample of clusters and groups of galaxies cov- 
ers the redshift range from z = 0. 1 to redshift z ~ I with a 
median redshift of z = 0.47. The median temperature of the 
clusters is ~ 2 keV, and the median Msoo mass 9 x lO'-' M© 
(based on luminosity scaling relations). With our ~ 10 ks 
XMM-Newton observations we are thus able to effectively 
probe the cluster/group transition regime practically at all 
redshifts up to z « 1 . 
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- We provide a preliminary, simplified calculation of the sur- 
vey sky coverage which does not require extensive Monte 
Carlo simulations. Using this calculation we characterize 
our cluster sample by its logA^ - log 5 relation. We find 
good ag reement with the re lations established by the RDCS 
surve y ( Rosati et"al] 19981). 400deg^ survey dBurenin et al.l 
[2007; 'Vikhlinin et all l2009l) and the XMM-LSS project 
(iPacaud et al. 2007^. 

- We carried out first comparisons with optical studies 
available from the Southern Cosmology Survey (SCS, 
iMenanteau et al.l l2009, 2010). In this preliminary investiga- 
tion we find the SCS photometric redshifts to be biased low 
by ~ 20% with respect to our estimates (both photometric 
and spectroscopic, where available). We find a discrepancy 
between the X-ray and optical mass estimates, with optical 
masses being significantly higher. We compare our masses to 
weak lensing mass measurements available for 13 clusters in 
our sample from Mclnnes et al.l (l2009i) . Although the weak 
lensing mass uncertainties are large, there is no statistical in- 
consistency between the two mass estimators. 

The presented results illustrate the potential of medium- 
deep. X-ray surveys to deliver cluster samples for cosmologi- 
cal modelling. These samples then in combination with avail- 
able multi-wavelength data (particularly in optical, near-infrared 
and SZE) will allow us to probe the dependence of the selection 
functions on relevant cluster observables and provide thus an im- 
portant input for upcoming large-area multi-wavelength cluster 
surveys. 
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Appendix A: Ancillary information 

A.1. Quality flags 

In this section we provide useful ancillary data for our clusters 
described using several X-ray quality flags and diagnostic pa- 
rameters compiled in Table lA.ll Here is the description of the 
table's columns: 

- ID: the cluster identification number 

- BCS field: the identification number of the BCS field, on 
which the cluster is lying. Some clusters lie on two or more 
tiles, in those cases we provide the name of the tile with the 
largest overlap region. 

- XMM OBSID: The official identification number of the 
XMM-Newton pointing containing the cluster If the cluster 
lies in two (or three) adjacent observations we provide the 
OBSID of the pointing which provides the best constraint on 
the cluster flux (typically the one where the cluster is at the 
smallest off'-axis angle). 

- Internal field ID: XMM-BCS internal field ID used in the 
text as shorthand for the OBSID. See also Table[T]for the full 
cross-listing of field IDs. 

- flagHc^ The hot chip flag is a four character string, with the 
characters being either T for "true" or F for "false". The sig- 
nificance of the characters: 

1 . character: Does the observation have a hot MOS2 CCD#5? 

2. character: Does flie cluster lie on the MOS2 CCD#5? 

3. character: Does the observation have a hot MOS 1 CCD#4? 

4. character: Does the cluster lie on the MOSl CCD#4? 
For the problematics of the hot chips see Sect. 13.1.11 

- SNRxflux: The flux estimation significance determined as 
^soo/cfsoo' where F500 is the source flux in the rgoo aper- 
ture, and (Tf^f^ is its error (including shot noise and 5% back- 
ground modelling uncertainty. Sect. 13.2b . 

- flaginst: The instrument flag equals if the physical param- 
eters of the source were obtained using the PN and both 
MOS cameras. If flaginst - 1 only PN could be used and 
if flaginst - 2 only the combination of the two MOS cameras 
was utilized. 

- Qpj^j: The automatic plateau fit quality flag for the PN 
growth curve s. Th e se fla gs have the same meaning as 
in Bohringer e t al.l (l2000l) . In summary, as described in 
Sect. 13.21 we fit a line to the growth curve between Tpi„, and 
the outer extraction radius. The flag describes the quality of 
this fit by calculating the ratio of the predicted count rate 
from the linear fit to the expectation, if the plateau was con- 
stant and equal to the estimated plateau flux. Q™^ = 1 : the 
growth curve shows neither significant increase nor decrease 
outside Tpiai- This value is assigned if the linear extrapola- 
tion does not differ by more than 0.8% per bin from the 
constant value. Q™^ - 2: marks a declining curve (decline 
> 0.8%/bin). A decline can occur if the background model 
(determined from a fit to the whole field) slightly overesti- 
mates the local background. In this case we attempt to es- 
timate the plateau level from the 3 bins closest to r;,/„,. If 
the final fit is acceptable (no significant residual decline), the 
plateau is accepted and assigned this quality flag. Q™^ = 3 

and Q™j - 4: in case that the plateau is rising an attempt 
is made to iteratively exclude the outermost bins and in a 
second step also the innermost bins. This procedure helps in 
correcting an outer rise of the growth curve due to a neigh- 
boring source and if necessary by skipping over a few bins 
if the curve fluctuates in the radial range close to r^/a,. If 



this procedure converges, after the exclusion of the outer- 
most bins the plateau is accepted and flagged with Q™^ = 3. 
If the procedure converges, but it required also the second 
step of excluding the innermost bins we assign Q™ - 4. If 

Qpiat ~ ^' plateau is rising and the increase could not be 
corrected for by the above described procedure. If there are 
only two or less radial bins outside the plateau radius can not 
be established and we assign Q™^ - 9. 

Qplat - ^ mark generally good quality plateaus (naturally, 
the lower the flag the better). Q™^ = 5 is a serious warning 
and Q™j = 9 is not recommended to be used at all. In fact, 
for the parameters in Table |2] we do not use plateaus with 
this flag with the exception for the systems ID 476 and 139 
where an alternative solution is not available due to their 
significant blending (Appendix lA.2b . 

- Q^OS; The same as Q™ but applied to the MOS1+MOS2 
growth curve. 

- Qg^: Visual flag set considering the overall quality of the PN 
growth curve solution (taking into account the presence of 
chip gaps, anomalous background, potential contamination 
etc.). Value equal to 1 is the best (no problems), equal to 3 the 
worse (to be considered as a warning). If the given detector 
was not included in the analysis (as marked by flaginst) Qg^ 
is set to 9. 

- Q^?^: The same as Q™ but appHed to the MOS1-HMOS2 
growth curve. 

- Qtot: This flag is a non-quantitative ancillary X-ray flag. 
It is assigned during a visual inspection of each source, but 
takes into consideration also all above flags. Sources with 
this flag equal to 1 (best) and 2 (only mild warnings) have 
high quality X-ray photometry measurements. Flag equal to 
3 signifies problems with the X-ray photometry e.g. due to 
problematic locations on the detectors, bad quality plateaus, 
source blending etc. In the present work, we have included 
also Qtot = 3 sources in all analyses since they constitute 
only about ~ 10% of the sample. 



Table A.l. Cluster ID and proper name cross-reference table. We also list the BCS field name for each cluster and additional ancillary X-ray quality flags (see text for 
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A.2. Notes on individual sources 

Some of the identified clusters required individual treatment and 
in this section we provide notes for these cases: 

- ID Oil: this high redshift cluster lies on the heavily flared 
field F03 with no quiescent period. Therefore, the field was 
not used for the sensitivity function calculation and the 
logA^ - log 5. The double component background model 
accounts in principle in the first approximation for the en- 
hanced background and therefore we provide the basic X-ray 
parameters for this cluster The diagnostic flags (Table lA.U 
indicate that the growth curve solution is quite reliable, but 
due to the flaring all physical parameters should be treated 
with caution. 

- ID 038: This source consists of two completely overlapping 
systems, one with photometric redshift z = 0.39 + 0.05 and 
the second with z = 0.74 + 0.07. Since there is no direct 
way to disentangle the contribution of the two sources, we 
will assume that all the flux comes from the more nearby 
system. In this case, the estimated physical parameters are 
upper limits. 

- ID 070: is a nearby cluster with large extent and measured 
flux. It lies on a hot MOS2 CCD#5 and due to its extent it is 
impossible to obtain a background area on this chip uncon- 
taminated by the source emission. Therefore we can not use 
the procedure described in Sect. 13.1.11 where we fit a dou- 
ble component model to the hot chip independently from the 
rest of the field. Instead we discard the data from this chip 
completely. Consistently, we use the MOS 1 response matrix 
instead of the MOS2 one for further analysis. 

- ID 081: in this distant cluster (photo-z= 0.85) the BCG is 
off'set by ~ 19" from the X-ray centroid. The emission from 
this system might be contaminated by an intervening AGN 
(the X-ray center is coincident with a bright galaxy which 
could harbor one). 

- ID 109: due to the limited depth of the available optical data 
we can provide only very tentative redshift estimate for this 
system. 

- ID 139 and 476: We detect two nearby, high significance 
extended sources in this region (~ 1.3 arcmin apart). The 
systems are confirmed as independent also in redshift space 
by our spectroscopic measurements (ID 476 at z = 0. 102 and 
ID 139 at z = 0.169). In order to measure the flux of each 
cluster we excise the other source. Due to their proximity, 
however, full deblending is not possible and therefore both 
fluxes are likely overestimated. The analysis of the sources 
is further complicated by the presence of a very bright X-ray 
point source at ~ 2 arcmin distance from the clusters and 
high quiescent soft proton contamination in the PN camera 
(39% background inc rease, Sect.|2ll . 

The cluster catalog of lBurenin et al.l (l2007h based on ROSAT 
data includes a source with a center roughly between the two 
systems (i.e. very likely misclassified as a single cluster due 
to the limited resolution of ROSAT). 

- ID 275: also lies on a hot M0S2 CCD#5. Its detection likeH- 
hood is completely dominated by the MOS2 detection, how- 
ever the source is not flagged as spurious based on the criteria 
described in Sect. 13.1.11 because it would be above the de- 
tection threshold even without the MOS2 data. In this case, 
the background modeling of the hot chip was possible and 
this background was used in the subsequent growth curve 
analysis. 



A.3. Galactic counterparts 

All NED galactic counterparts within 16" from the X-ray cen- 
troids of our clusters are summarized in Table IA.2I (see also 
Sect. El. 

A.4. Wavelet detections 

We have also carried out a detection using a wavelet de- 
tection algorithm deve l oped for the COSMOS project by 
iFinoguenov et al.l (l2007l l2009l) . Every cluster presented in our 
current sample has also been confirmed by this approach. In 
addition, we have identified five systems (Table IA.3I ) with a 
wavelet detection but no SAS-based detection in any setup. All 
five systems are coincident with significant galaxy overdensities. 
We find many interloping X-ray point-sources in these systems 
(a potential source of misclassification in the SAS detections). 
Even after conservative point source removal, residual contam- 
ination makes all the estimated X-ray parameters for these sys- 
tems highly uncertain. These detections are not included in the 
statistical description of the sample (e.g. the logA^ - log 5 rela- 
tion) and are listed here only for completeness (in Table lA.3T l. 
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Table A.2. Galaxies identified in the NED database to be within 16" from the X-ray center. We list spectroscopic redshifts where 
available - in both cases the identified galaxies are BCGs. Redshift reference: " I Jones et al ] (12004 . 



iU 


Object Name 


K.A. {aeg) 




redshift 


separation 


HQ 4 


ArMUKs(r>J) r>/oz/jU.lU-5DUU12.1 


Q iiC /I A 


-J J. 11 /4 




1 A" 

1.9 


Ujy 


2MAJ5A J/jiyi / 12-3519254 




C C QO /I C 

-55.3243 




A 

U.5 




2JvlAoA J2jlyU212-jj2jl9D 


349. / joo 


-55.3085 




1.3 


070 


ZMAbX J2j22j092-5416086 


350.6289 


-54.2691 




0.8" 


094 


APMUKS(BJ) B232918.62-552918.2 


353.0196 


-55.2122 




2.5" 


127 


2MASX J23272468-5503589 


351.8528 


-55.0664 




9.6" 


150 


2MASX J23300047-5437069 


352.5019 


-54.6187 


0.177" 


1.5" 


152 


2MASX J23294006-5447220 


352.4168 


-54.7895 




3.1" 


227 


APMUKS(BJ) B23 1920.30-554137.6 


350.5444 


-55.4194 




4.4" 


268 


APMUKS(BJ) B232326. 14-554657.3 


351.5618 


-55.5074 




14.0" 


476 


2MASX J232540 15-5444308 


351.4173 


-54.7419 


0.101" 


3.1" 


511 


APMUKS(BJ) B232929.68-545847.0 


353.0645 


-54.7035 




10.8" 


547 


2MASX J23241957-5525494 


351.0816 


-55.4303 




0.6" 



Table A.3. Physical parameters for the lower quality detections. 



ID 


R.A. (J2000) 


Dec (12000) 


z 




^ phiT 


^5m 


^500 


^500 


7^500 




^^500 


M200 




(deg) 


(deg) 


(photo.) 


(ai-cmin/r^fli,) 


(Ur'"' ergs"' cnr-) 


(kpc) 


(10""ergs"'cm"^) 


(10«ergs-') 


(keV) 


(W" Mo) 


(lO" Mo keV) 


(10'-' Mo) 


355 


352.8205 


-54.7064 


0.31 ±0.05 


2.2/1.1 


3.20 ± 0.94 


541 ±55 


2.66 ± 0.79 


0.8 ± 0.2 


1.6 ±0.5 


6.2 ± 1.9 


0.6 ± 0.4 


8.5 ±2.6 


534 


351.1993 


-55.4128 


0.31 +0.02 


1.4/0.7 


2.18± l.ll 


523 ± 64 


2.22± 1.13 


0.7 ± 0.4 


1.5 + 0.5 


5.6 ± 2.0 


0.5 ± 0.4 


7.6 ±2.8 


536 


354.2994 


-55.3514 


0.69 ± 0.07 


0.8/0.7 


0.83 ± 0.23 


487 ± 49 


0.85 ± 0.24 


1.8 ±0.5 


1.9 ±0.6 


7.1 ±2.1 


0.8 ±0.5 


10.0 ±3.0 


540 


350.1461 


-54.4628 


0.24 ± 0.02 


1.3/0.6 


1.66 ±0.46 


474 ± 48 


1.75 ±0.48 


0.3 ±0.1 


1.2 ±0.3 


3.8 ± 1.2 


0.2 ±0.1 


5.2 ± 1.6 


541 


350.5048 


-54.3296 


0.59 ± 0.06 


0.9/0.7 


1.06 ±0.30 


506 ±51 


1.08 ± 0.30 


1.5 ± 0.4 


1.8 ± 0.5 


7.0 ±2.1 


0.8 ± 0.5 


9.9 ± 3.0 
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Appendix B: Test of the sensitivity function from 
preliminary l\/lonte Carlo simulations 

The determination of the survey selection function is a crucial 
requirement for the cosmological modeling of the cluster sam- 
ple, scaling relation studies etc. Due to the complex nature of 
the extended source detection, this question can be properly ad- 
dressed only by detailed Monte Carlo simulations. In the present 
work we utilized a simplifying approach that allowed us to get 
a first estimate of the sensitivity functions and the recovered 
log N-logS relation (Sect. |52]and|l3]l- 

T he software for the Monte Carlo simulations (Miihleggerl 
l2oToh is in an advanced development stage which allows us to 
carry out a preliminary test of our simplified approach. 

The simulation pipeline uses the survey fields themselves 
and injects mock beta model clusters into the observations at 
random positions across the field-of-view. The field is then pro- 
cessed with the detection pipeline. The process is repeated on 
a grid of cluster fluxes and core radii and the cluster detection 
probability is derived as a function of these parameters. The use 
of real observations instead of model backgrounds allows us to 
derive a realistic selec tion function. Th e simulation software is 
described in detail in MuhleggeddloTol) . 

Simulations are currently available fo r a subset of the XMM- 
Newton Distant Cluster Project (XDCP. Bohringer et ai]l2005t 
[passbender 2008; Fassbender et al. 2011a, submitted) fields. 
From these fields we selected 3 observations (XMM OBSIDs 
0104860201, 0111970101, 0112551101) which have similar 
depth to our survey fields (e.g. cleaned exposure times ~ 10 ks 
and enough area unaff'ected by the central source to safely as- 
sess the background). We processed these fields with our de- 
tection pipeline and calculated the point source sensitivity func- 
tion and the scaled extended sensitivity function as described in 
Sect. 15.21 The comparison with the sensitivity function derived 
from the simulations are displayed in Fig. IB. II The simple cal- 
culation (black curves) akeady matches the realistic calculation 
(red curves) very well, capturing also the transition parts of the 
curve. The curves from simulations include the effect of incom- 
pleteness of the output catalogs. The red curves in Fig. IB. II are 
calculated for a 50% completeness level (c=0.5). The complete- 
ness of our cluster catalog can be assessed only by simulations, 
but is certainly higher than 50%. This means that the prelimi- 
nary analytic sky coverage function overestimates the sky cov- 
erage. The use of the true sky-coverage function would lead to 
an increase of the weighting factor in Eq.|5]and would move the 
points in Fig. [12] in the relevant flux range slightly higher The 
sensitivity function for a 90% completeness scenario is plotted 
in green and as expected yields a much smaller area for the given 
flux. 

Additional subtle effects slightly influence this comparison, 
e.g. leading to different normalizations of the two curves in the 
saturated high-end part: 1) All the fields have a bright source 
in the center of the field-of-view, which has to be excised. The 
excision is treated slightly differently in the simulations and in 
the simple calculations leading to slightly different total geomet- 
ric areas. 2) The simulated curves were calculated for the single 
band detection scheme in the 0.35 - 2.4 keV while our analytic 
solution is for the 0.5 - 2 keV band. The fluxes were converted 
to the 0.5 - 2 keV band, but detection in these different energy 
ranges could cause slightly different completeness and contami- 
nation fractions. 

We conclude, that our first-order approach yields a good 
description of the sensitivity function for a 50% completeness 
level. To estimate the completeness of our sample Monte Carlo 




10"" 

t;,(0.5-2.0 keV) [erg s-'cm-=] 



Fig.B.l. Sky-coverage for extended sources in three XDCP 
fields from Monte Carlo simulations at the 50% (red curve) and 
90% (green) completeness level. The black solid curve shows 
the sky-coverage calculated by scaling the point source curve 
(dashed) with an offset factor of 2.4 (see Sect. 15.2b . The simple 
scaling is shown to be a good first order description of the ex- 
tended source sensitivity function. See AppendixiBjand Sect. 15.21 
for details. 



simulations are needed. The sensitivity functions from Sec 
provides sufficient precision for present applications and the pre- 
liminary log - log S is already in good agreement with previ- 
ous findings. The described simulation pipeline will be applied 
to the whole XMM-BCS survey in subsequent work and the re- 
alistic selection function will be utilized for further analysis and 
modelling of the final cluster sample. 



Appendix C: Comparison with the XMM-LSS survey 

The first part of the XMM-LSS survey (the initial 5 deg^, 
Pier re et al. I l2007t iPacaud et al.ll2006l l2007h offers an excellent 
match to our survey not only with respect to the area, but also 
to the typical depth (having only slightly higher average expo- 
sure times). Since the XMM-LSS project has already carried out 
Monte Carlo simulations to calibrate its detection and source- 
characterization pipeline, we make here an effort to compare 
results derived from our XMM-BCS pipeline with results pub- 
lished by XMM-LSS. 



C.1. Cluster detection comparison 

A full comparison of the source detection pipelines would be 
only of limited use and is cmTently impossible since only a 
small part of the XMM-LSS extended sources have been thus 
far spectroscopical ly confirmed (the so-called CI sample of 
Pacaud et al. 2007) Pi Therefore, we restrict ourselves to the re- 
analysis of the CI sample. 

We downloaded all the XMM-LSS fields with CI detec- 
tion£3 and fully reanalyzed them with the XMM-BCS pipeline. 



Catalog available at : 
heasarc . gsf c . nasa . gov/W3Browse/ all/xitmlssoid . html 

'5 XMM OBSIDs: 0037980301, 0037980701, 0037981001, 
0037981101, 0037981201, 0037981501, 0037981601, 0037981801, 
0037982501, 0037982601, 0109520201, 0109520301, 0109520601, 
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We confidently detected all the CI clusters and they are among 
our highest ranked extended source detections. 

In Fig. IC.ll we compare their detection and extent 
likelihoods with their respective XMM-LSS variants 
(SB_Detect_Likelihood and SB_Extent_Likelihood). 
Both sets of parameters exhibit a strong correlation, showing a 
good consistency between both detection approaches (XMM- 
LSS uses a single band wavelet detection scheme). The scatter 
between the parameters is caused by differences in the data 
reduction process, background estimation and source detection 
algorithms. 

The CI sample is defined by SB_Detect_Likelihood> 32, 
SB_Extent_Likelihood> 33. We fit a hnear relation in the 
two log-log planes and use these cuts to convert the XMM-LSS 
thresholds to our parameters obtaining: det_ml> 16.4 (equiva- 
lent to ~ SAcr detection in our scheme) and ext_ml> 8.3 (i.e. 
~ 3. 7 cr extent significance). 

C.2. X-ray photometry comparison 



In Fig. IC.2l we compare the fluxes in the 0.5 - 2 keV band and 
0.5 Mpc aperture measured by the XMM-LSS and by us using 
the growth curve method (Sect. [J!2T i. Being interested in the flux 
estimation we have fixed the cluster redshifts to its spectroscopic 
value provided by XMM-LSS. We choose not to use the infor- 
mation on the spectroscopic temperature, but rather we estimate 
it from the scaling relation as we do for the XMM-BCS sample 
(see SectlOi. 

The 0.5 - 2 keV fluxes in a 0.5 Mpc aperture are compared 
in Fig. IC.2I where the residuals in the bottom panel are defined 

as Afx = f^MM-LSS _ j:XMM-BCS -^g f^^^ jj^g fl^j^-g^ jjj ^gj.gg. 

ment with the mean relative difference of ~ 0.6% and a scatter of 
~ 17% (averaged over the whole sample). If we split the sample 
by the median flux (fx w 3.9 x 10"'"^ erg cm"^ s"') the XMM- 
LSS fluxes are on average by ~ 4% higher on the faint end while 
being slightly lower by roughly the same amount on the brighter 
end. A systematic difference of this magnitude would however 
cause only a ~ 1% shift in the temperature estimates and ~ 2% 
shift on mass. An analogous comparison for the bolometric lu- 
minosity in the rjoo aperture is displayed in Fig. IC.3I Here the 
mean relative difference is ~ 6.6% and a scatter of ~ 19% (the 
corresponding temperature/mass difference are 2%/3%). 

The overall agreement is encouraging, if we take into ac- 
count that the two pipelines utilize principally different ap- 
proaches to the flux measurement. XMM-LSS utilizes a beta 
model fit to the cluster's surface brightness integrated out to a 
fiducial radius, while our method is completely non-parametric 
(except for a typically small extrapolation factor if the required 
aperture is larger than the range where the cluster emission is 
detected directly). Background estimation in both approaches is 
also markedly different. 

Since we did not use the information on the spectroscopic X- 
ray temperature, it is interesting to note, that the mean tempera- 
ture residuals (spectroscopic compared to our estimates from the 
scaling relations) are < 2% with a standard deviation of ~ 23% 
(i.e. comparable to measurement errors). Although the error bars 
are large, this agreement indicates that the L-T scaling relation 
and its evolution adopted in this work (iPratt et al.ll2009h are suit- 
able for cluster samples drawn from surveys of this type. 

The cluster mass is not a direct observable in either of the 
two surveys. XMM-LSS gives rough estimates based on their 



spect roscop i c mea surement and beta model fit using the relation 
from Ettori (2000). Our estimates, using the L - M relation of 
Pratt et al. (2009), give on avera ge a lmost 40% higher masses. 
We note however, that the lEttoril (12000 ) relation was derived be- 
fore the advent of XMM-Newton and Chandra. It is thus not 
well representing our current understanding of cluster mass esti- 
mation. 

Finally, we also check the consistency of the beta model fits 
between the two pipelines. Since the core radius ^core and the /3 
exponent of the beta model are strongly degenerate, especially 
for the case of low counts profiles, our fitting procedure keeps 
/3 fixed to the canonical value of 2/3. The XMM-LSS pipeline 
carries out fits with both rcore and /3 as free parameters. Despite 
this difference, we find good agreement between the estimated 
core radii (Fig. lC.4| ). 




XMM-BCS fx [lO '" erg s ' cm '^] 

Fig. C.2. Comparison of measured X-ray fluxes for the CI sub- 
sample of the XMM-LSS survey in the 0.5 - 2 keV band and a 
0.5 Mpc aperture (iPacaud et al.l2007t y-axis) and the fluxes mea- 
sured by our pipeline (x-axis). The red line marks equality. The 



bottom panel shows the residuals Afx = f^ 
See Sect. lC.2l for details. 



■XMM-BCS 



0111110301, 0111110401, 0112680101, 0112680201, 0112680301, 
0112680401, 0112680501, 0147110101, 0147110201. 
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Fig. C.l. Compa rison of detection ( left panel) and extent likelihoods (right panel) between our pipeline (x-axis) and the XMM-LSS 
pipeline (y-axis. lPacaud et al]|2007h . The derived likelihoods are well correlated and the red line shows the best fit relations. 
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Fig. C.3. Comparison of bolometric X-ray luminosities for the 
CI subsample of t he XMM-LSS survey in the rsoo aperture 
jPacaud et alj l2007l y-axis) and the luminosities measured by 
our pipeline (x-axis). The red line marks equality. The bottom 
panel shows the residuals ALx = L^^m-lss _ ^xmm-bcs 
Sect. IC.2lfor details. 



Fig. C.4. Beta model core radii for the XMM-LSS CI sample as 
estimated by our pipeline (x-axis) and by the XMM-LSS. The 
red line marks equality. The core radii are typically highly un- 
certain given the relatively low photon statistics. Despite this, 
the agreement between the two estimates is good. Note that the 
XMM-LSS values are fitted with the beta value as a free param- 
eter, while we fix its value to 2/3. 



